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Abstract. 

We present the hrst measurements of the Probability Distribution Function (PDF) of galaxy fluctuations in the 
four-passes, first-epoch VIMOS- VLT Deep Survey (VVDS) cone, covering 0.4x0.4 deg between 0.4< z<1.5. We 
, show that the PDF of density contrasts of the VVDS galaxies is an unbiased tracer of the underlying parent 

distribution up to redshift 2=1.5, on scales R=8 and 10 /i _1 Mpc . The second moment of the PDF, i.e., the 
rms fluctuations of the galaxy density field, is with good approximation constant over the full redshift baseline 
investigated: we find that, in redshift space, as for galaxies brighter than M C B = —20 + 51og/i has a mean value 
of 0.94±0.07 in the redshift interval 0.7<z<1.5. The third moment, i.e., the skewness, increases with cosmic 
time: we find that the probability of having underdense regions is greater at 2~0.7 than it was at 2~1.5. By 
comparing the PDF of galaxy density contrasts with the theoretically predicted PDF of mass fluctuations we 
infer the redshift-, density-, and scale-dependence of the biasing function b(z,8,R) between galaxy and matter 
overdensities up to redshift z = 1.5. Our results can be summarized as follows: i) the galaxy bias is an increasing 
function of redshift: evolution is marginal up to z~0.8 and more pronounced for z^0.8; ii) the formation of bright 
galaxies is inhibited below a characteristic mass-overdensity threshold whose amplitude increases with redshift 
and luminosity; iii) the biasing function is non linear in all the redshift bins investigated with non- linear effects 
of the order of a few to ~ 10% on scales > 5 /i _1 Mpc. By subdividing the sample according to galaxy luminosity 
and colors, we also show that: iv) brighter galaxies are more strongly biased than less luminous ones at every 
redshift and the dependence of biasing on luminosity at z~0.8 is in good agreement with what is observed in the 
local Universe; v) red objects are systematically more biased than blue objects at all cosmic epochs investigated, 
but the relative bias between red and blue objects is constant as a function of redshift in the interval 0.7< z<1.5, 
and its value (b rel ~ 1.4) is similar to what is found at 2~0. 



Key words, cosmology: deep redshift surveys — cosmology: theory — large scale structure of the Universe — galaxies: 
distances and redshifts — galaxies: evolution — galaxies: statistics 



1. Introduction 

The understanding of how matter structures grow via 
gravitational instability in an expanding Universe is quite 
well developed and has led to a succe s sful and predictive 
theor etical framework (e.g., IPeeblesL Il98f)t iDavis et all 

One of the most critical problems, however, is to un- 
derstand the complex mechanisms which, on various cos- 
mological scales, regulate the formation and the evolution 
of luminous structures within the underlying dark-matter 
distribution. Its solution ultimately relics on the compre- 
hension of the "microscopic" physics which describes how 
the baryons fall, heat-up, virialize, cool and form stars in 
the potential wells generated by the dominant mass com- 
ponent of the Universe, i . e., the non-baryonic dark mat- 
ter (e.g., IWhite fc R.eesL Il978h . A zero-order, minimal 
approach to this investigation consists in characterizing 
"macroscopically" the cosmological matter fluctuations in 
terms of a reduced set of fundamental quantities, essen- 
tially their positions and mass scales, and in studying how 
the respective spatial clustering and density amplitudes 
relate to the corresponding statistics computed for light 
fluctuations. This comparison sche me is generally referre d 
to as matter-galaxy biasing (e.g., iDekel fc Lahavl Il999) . 

An operational definition of bias is conventionally 
given in terms of continuous density fields by assuming 
that the local density fluctuation pattern traced by galax- 
ies (S g ) and mass (S) are deterministically related via the 
"linear biasing scheme" 



Sg(6) = bS 



(1) 



where the co nstant "slope" b is the biasing parameter 
feaiseilll984L 

This specific formulation, however, represents a very 
crude approximation which is not based on any theoret- 
ical or physical motivation. It is obvious, for example, 
that such a model cannot satisfy the physical requirement 
S g (— 1) = —1 for any arbitrary value b ^ 1. In particular 
the b iasing process cou ld be non local (e.g.. Iciatelan et all 
Il998l) . sto chastic (e.g., [ Dekel fc LahavMl999() and non lin- 
ear (e.g., iMo fc White! . 199fift . Moreover, both theory and 
numerical simulations predict that the bias grows mono- 
tonic ally from the p r esent cosmic epoch to high redshifts 
(e.g.. IDekel fc Ree.4 li^ If71 Il99l iMo fc Whitfl.ll993 

^ii|2QQ2). 



iTegmark fc Peebles! ll998TlBasilakos fc Plionisl 12001 

From a theoretical perspective, light does not fol- 
low the matter distribution on sub-galactic scales, where 
nearly 90% of dense, low-mass dark matter fluctuations 
(M~ 1O 7 -1O 8 M 0) failed to form stars and to become 
galaxies (e.g.. iKlvnlin et, all [Oil iMoore et all [jjil 
iDalal fc Kochanekt l2002|) . A difference in the spatial dis- 
tribution of visible and dark matter is predicted also 
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on galactic scales, since the radial scaling of density 
profil es of dark matter halos l|Navarro. Frenk fc WhiteL 
Il997h differs from the three-dimensional radial distribu- 
tions of light (Sersic and Freeman laws). Galaxy biasing 
is theoretically expected also on cosmological scales. In 
particular, simulations of the large-scale structure pre- 
dict the existence of a difference in th e relative distri- 
butio n of mass and dark halos (e.g., ICen fc Ostrikerl 
Baglal Il998t iKravtsov fc Klypirl Il999h or galaxies 
Evrard. Summers. & Oavis, I 9941: iBlanton et all 



EM 

(e.g., 

l2000t IKavo. Taruva fc Sutol l200l|) . Various physical 
mechanisms for biasing have been prop osed, such as , 
for example, the p eaks-biasing scheme l)Kaiserl Il984l 
iBardeen et all [l986^ . the probabilistic biasing approach 
dColesL I199.C or the biasing models derived in the con- 
text of the extended Press fc Schechter a pproximation 
llMo fc WhiteL ll99fitlMatarrese et allll998h . 

Turning to the observational side, the fact that, in 
the local Universe, galaxie s cluster differently a ccord- 
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can simultaneously trace the underlying distribution of 
mass, and that galaxy biasing not only exists, but might 
also be sensitive to various physical processes. As a mat- 
ter of fact, redshift information which recently became 
available for large samples of galaxies have significantly 
contributed to better shaping our current understanding 
of galaxy biasing, at least in t he local Universe. T he anal- 
ysis of th e power spectrum 



ill . I2002I) and bi- 

spectrum JVerde et alll2002t) o f the 2 dF Galaxy Redshift 
Survey (2dFGRS IColless et all l200l|) consistently shows 
that a flux-limited sample of local galaxies ( z<0.25), opti- 
cally selected in the 6j-band [bj < 19.4), traces the mass, 
i.e., it is unbiased, on scales 5<R(/i~ 1 Mpc )<30. 

The galaxy correlation fun ction has been measur ed up 
to redshift ~1 by the CFRS llLe Feyre et all Il99^ . and 
by the CNOC llOarlberg et aim OOO) surveys giving con- 
flicting evidence on cluster ing amplitude and bias evolu- 

More recently, the analysis of 



tion (see 



the first season PEEP2 data l|Ooil et all Eool seems to 



indicate that a combined R-band plus color selected sam- 
ple is unbiased at z~l. On the contrary, measurements of 
the clustering JSteidel et allll998tlGiavalisco et allll998l 
iFoucaud et alll200. c jl) or of the ampli t ude o f the count- 
in-cell fluctuations l)Aldeberger et all Il998l) of Lyman- 
break galaxies (LBGs) at z~3 suggest that these objects 
are more highly biased tracers of the mass density field 
than are galaxies today. Higher redshift domains have 
been probed b y usin g photometric redshift information 
l|Arnouts et all Il999l). or compilat ion of heterogeneous 
samples l(Magliocchetti et all l2000l) . Again, the cluster- 
ing signal appears to come from objects which are highly 
biased with respect to the underlying distribution of mass. 

While there is general observational consensus on the 
broad picture, i.e., that biasing must decrease with cosmic 
time, the elucidation of the finer details of this evolution 
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as well as any meaningful comparison with specific the- 
oretical predictions is still far from being secured. Since 
clustering depends on morphology, color and luminosity, 
and since most high redshift samples have been selected 
according to different colors or luminosity criteria, it is 
not clear, for example, how the very different classes of 
objects (Ly-break galaxies, extremely red objects or ul- 
traluminous galaxies), which populate different redshift 
intervals, can be considered a uniform set of mass trac- 
ers across different cosmic epochs. Furthermore, one must 
note that the biasing relation is likely to be nontrivial, 
i.e., non-lin ear and scale dependent , especially at high red- 
shift (e.g., ISomerville et all 1200 ill 

Only large redshift surveys defined in terms of uni- 
form selection criteria and sampling typical galaxies (or 
their progenitors), rather than small subclasses of pecu- 
liar objects, promise to yield a more coherent picture of 
biasing evolution. In particular, the 3D spatial informa- 
tion provided by the V IMOS-VLT Deep Survey (VVDS, 
iLe Fevre et all l2005al hereafter Paper I) should allow us 
to investigate the mass and scale dependence, as well as to 
explore the time evolution of the biasing relation between 
dark matter and galaxies for a homogeneous, flux-limited 
(I<24) sample of optically selected galaxies. 

The intent of this paper is to provide a measure, on 
some characteristic scales R, over the continuous redshift 
interval 0.4< z<1.5, of the local, non-linear, deterministic 
biasing function 

b = b{z,6,R). (2) 

The goal is to provide an observational benchmark to 
theories predicting the efficiency of structure formation, or 
semi-analytical simulations of galaxy evolution. The prob- 
lem, however, is to find an optimal strategy to evaluate 
the biasing function in the quasi-pencil-beam geometry 
of the first-epoch VVDS survey. At present, the angular 
size of the first-epoch VVDS redshift cone (~ 0.5deg 2 ) 
does not allow us to constrain the biasing function us- 
ing high order moments of the galaxy distribution (for 
example the 3-point correlation function). Moreover, we 
cannot determine the biasing function simply regressing 
the galaxy fluctuations (5 g ) versus mass fluctuations (5). 
The fundamental limitation preventing such an intuitive 
comparison is evident: it is easy to "pixelize" the survey 
volume and to measure the galaxy fluctuations in each 
survey cell, but, since the VVDS is not a matter survey, it 
is much less straightforward to assign to each cell a mass 
density value. Thus, in this paper, we take an orthogonal 
approach and we infer the biasing relation S g = 5 g (S) be- 
tween mass and galaxy overdensities from their respective 
probability distribution functions (PDFs) f(S) and g{S g ): 
assuming a one-to-one mapping between mass and galaxy 
overdensity fields, conservation of probability implies 

dS g (S) _ f(5) 
dS g(S g ) ' 

The advantage over other methods is that we can 
explore the functional form of the relationship S g = 



b{z, S, R)S over a wide range in mass density contrasts, 
redshift intervals and smoothing scales R without specify- 
ing any a-priori parametric functional form for the biasing 
function. 

In pursuing our approach, we assume that the PDF 
of matter overdensities f(6) is satisfactorily described by 
theory and N-body simulations. What we will try to assess 
explicitly, is the degree at which the measured PDF of 
the VVDS overdensities g(6 g ) reproduces the PDF of the 
underlying parent population of galaxies. The large size 
and high redshift sampling rate of the VVDS spectroscopic 
sample, together with the multi-color information in the 
B, V, R, I filters of the parent photometric catalog and the 
relatively simple selection functions of the survey, allow us 
to check for the presence of observational systematics in 
the data. In principle, this analysis helps us to constrain 
the range of the parameter space where first-epoch VVDS 
data can be analyzed in a statistically unbiased way and 
results can be meaningfully interpreted. 

The outline of the paper is the following: in §2 we 
briefly describe the first-epoch VVDS data sample. In 
§3 we introduce the technique applied for reconstruct- 
ing the three-dimensional density field traced by VVDS 
galaxies, providing details about corrections for various 
selection effects. In §4 we outline the construction of the 
PDF of galaxy overdensities and test its statistical rep- 
resentativity. We then derive the PDF of VVDS density 
contrasts and analyze its statistical moments. In §5 we re- 
view the theoretical properties of the analogous statistics 
for mass fluctuations. Particular emphasis is given to the 
problem of projecting the mass PDF derived in real space 
into redshift-perturbed comoving coordinates in the high 
redshift Universe. The method for computing the biasing 
function is introduced and tested against possible system- 
atics in §6. VVDS results are presented and discussed in 
§7. and compared to theoretical models of biasing evolu- 
tion in §8. Conclusions are drawn in §9. 

The coherent cosmological picture emerging from in- 
dependent observations and analysis motivate us to frame 
all the results presented in this paper in the context of a 
ACDM cosmological model with OT = 0.3 and — 0.7. 
Throughout, the Hubble constant is parameterized via 
h = i/o/100. All magnitudes in this paper are in the AB 
system (Oke & Gunn 1983), and from now on we will drop 
the suffix AB. 

2. The First-Epoch VVDS Redshift Sample 

The primary observational goal of the VIMOS-VLT 
Redshift Survey as well as the survey strategy and first- 
epoch observations in the VVDS-0226-04 field (from now 
on simply VVDS-02h) are presented in Paper I. 

Here it is enough to stress that, in order to mini- 
mize selection biases, the VVDS survey in the VVDS- 
02h field, has been conceived as a purely flux-limited 
(17.5 < I < 24) survey, i.e., no target pre-selection ac- 
cording to colors or compactness is implemented. Stars 
and QSOs have been a-posteriori removed from the final 
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redshift sample. Photometric data in this field are com- 
plete and free from surface brig htness selection effects, u p 
to the limiting mag nitude 1=24 l|Mc Cracken et al.ll2003ft . 

First-epoch spectroscopic observations in the VVDS- 
02h field were carried out using the VIMOS multi-object 
spectrograph l)Le Fevre et all l2003[) during two runs be- 
tween October and December 2002 (see Paper I). VIMOS 
observations have been performed using 1 arcsecond wide 
slits and the LRRed grism which covers the spectral range 
5500 < A (A) < 9400 with an effective spectral resolu- 
tion R ~ 227 at A = 75001. The accuracy in redshift 
measurements is ~275 km/s. Details on o bservations and 
data r eduction are given in Paper I, and in lLe Fevre et all 
(l2004h . 

The first-epoch VVDS-02h data sample extends over a 
sky area of 0.7x0.7 deg (which was targeted according to 
a 1, 2 or 4 passes strategy, i.e., giving to any single galaxy 
in the field 1, 2 or 4 chances to be targeted by VIMOS 
masks (see fig. 12 of paper I)) and has a median depth of 
about z~0.76. It contains 6582 galaxies with secure red- 
shifts (i.e., redshift determined with a quality flag>2 (see 
Paper I)) and probes a comoving volume (up to 2=1.5) of 
nearly 1.5 • 10 6 /i _3 Mpc 3 in a standard ACDM cosmology. 
This volume has transversal dimensions ~ 37x37 ft.~ 1 Mpc 
at z=1.5 and extends over 3060 /i _1 Mpc in radial direc- 
tion. 

For this study we define a sub-sample (VVDS-02h-4) 
with galaxies having redshift z<1.5 and selected in a con- 
tinuous sky region of 0.4x0.4 deg which has been homo- 
geneously targeted four times by VIMOS slitmasks. Even 
if we measure redshifts up to z~5 and in a wider area, 
the conservative angular and redshift limits bracket the 
range where we can sample in a denser way the underlying 
galaxy distribution and, thus, minimize biases in the re- 
construction of the density field (see the analysis in §4.1). 
The VVDS-02h-4 subsample contains 3448 galaxies with 
secure redshift (3001 with 0.4 < z < 1.5) and probes one- 
third of the total VVDS-02h volume. This is the main 
sample used in this study. 



3. The Density Field Reconstruction Scheme 

The first ingredient we need in order to derive the biasing 
relation 



S g = b{z,S,R)S 



(4) 



is a transformation scheme for diluting an intrinsic point- 
like process, such as the galaxy distribution in a redshift 
survey, into a continuous 3D ove rdensity field (see the re- 
view bv lStrauss fc WillickL fl995h . We write the dimension- 
less density contrast at the comoving position r, smoothed 
over a typical dimension R as 



5 g (r,R) 



Pg{r,R)- Pg 

~P~9 



(5) 



and we define (e.g., Irludsonl |l993) the smoothed number 
density of galaxies above the absolute magnitude thresh- 



old M. c as the convolution between Dirac's delta functions 
and some arbitrary filter 



S D (r 

p g (r,R,<M c ) = J2 — 



n) * 



|r — r 



S(ri,M c )$ z (m) 



(0) 



where the sum is taken over all the galaxies in the sample, 
S(r, M c ) is the distance-dependent selection (or incom- 
pleteness) function of the sample (see §3.2), ^ z {m) is the 
redshift sampling function (see §3.3) and F(\r\/R) is a 
smoothing kernel of width R. In this paper, the smooth- 
ing window F is modeled in terms of a normalized Top-Hat 
(TH) filter 



R 



(AnR 3 ) 



9 



(7) 



where O is the Heaviside function, defined as Q(x) — 1 for 
< x < 1, and <d(x) = elsewhere. 

Within this weighting scheme, shot-noise errors are 
evaluated by computing the variance of the galaxy field 



<r) = - 

Pg 



E 



F 



|r-ri| ' 



S{r u M c )®Mi) 



1/2 



(8) 



Note that all coordinates are comoving and that the 
mean density ~p~g~ depends on cosmic time. Since we observe 
an evolution of nearly a factor of two in the mean density 
of galaxies brighter than Ai* B (z = 0) from redshift z=0 
up to z=l (see Ilbert et al. (2005), hereafter Paper II), we 
compute the value of p~g~ at position r (corresponding to 
some look-back time t) with eq. © by simply averaging 
the galaxy distribution in survey slices r±R s where R s = 
400 h~ 1 Mpc. We note that conclusions drawn in this paper 
depend very weakly on the choice of Rs in the interval 
200< R s { h^Mpc )<600. 

In this paper, the density field is evaluated at positions 
r in the VVDS-02h volume that can be either random or 
regularly displaced on a 3D grid (see section 3.4). Even if 
we correct for the different sampling rate in the VVDS- 
02h field (by weighting each galaxy by the inverse of red- 
shift sampling function (<i> 2 (m))) we always select, for the 
purposes of our analysis, only the density fluctuations re- 
covered in spheres having at least 70% of their volume 
in the denser 4-passes volume. This in order to minimize 
the Poissonian noise due to the sparser redshift sampling 
outside the VVDS-02h-4 field. 

We also consider only volumes above the redshift 
threshold z t where the transversal dimension L of the 
first-epoch VVDS-02h cone is L(zt) >2R. As an example, 
for TH windows of size R=5(10) /i _1 Mpc we have Zt ~ 
0.4(0.7). Within the redshift range 0.4< z<1.5 the VVDS- 
02h field contains 5252 galaxies (of which 3001 are in the 
four passes region.) 

Note that we have characterized the galaxy-fluctuation 
field in terms of the number density contrast, instead of 
the luminosity density contrast, because the former quan- 
tity is expected to show a time-dependent variation which 
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is more sensitive to the galaxy evolution history (forma- 
tion and merger rates, for example). Moreover, as de- 
scribed in section §3.4, a robust description of the density 
field and a reliable determination of the PDF shape can 
be obtained only minimizing the shot noise component of 
the scatter; this is more easily done by considering galaxy 
number densities rather than galaxy light densities. 

The most critical elements of the smoothing process 
are directly readable in eq. 10: we must first evaluate 
galaxy absolute magnitudes at each redshift in the most 
reliable way, then specify the selection function and the 
redshift sampling rate of the VVDS survey. In the next 
sections we will describe how these quantities have been 
evaluated. 

3.1. The K-correction 

The absolute magnitude is defined as: 



AT 



5\ogd L (z,n) -K(z,SED) 



(9) 



where the suffixes r and o designate respectively the rest- 
frame band in which the absolute magnitude is computed 
and the band where the apparent magnitude is observa- 
tionally measured, and dh is the luminosity distance eval- 
uated in a given a-priori cosmology (i.e., using an appro- 
priate set of cosmological parameters l~i = [f2 m , ^a]). 

The correction factor K, which depends on redshift 
and the spectral energy distribution (SED), accounts for 
the fact that the system response in the observed frame 
corresponds to a narrower, bluer rest-frame passband, de- 
pending on the redshift of the observed object. A complete 
description of the application of this transformation tech- 
nique to VVDS galaxies is detailed in Paper II. 

The estimate of the galaxy absolute luminosity is thus 
affected by the uncertainties introduced by probing red- 
shift regimes where the k — correction term cannot be 
neglected. Anyway, using mock catalogs simulating the 
VVDS survey, we have shown (see Fig. A.f in Paper II) 
that the errors in the recovered absolute magnitude are 
significantly smaller in the B band (erg = 0.08) than in 
the I band (crj = 0.17). Thus, in what follows we will 
use absolute luminosities determined in the B-band rest- 
frame. 



3.2. The Radial Selection Function 

Since our sample is limited at bright and faint apparent 
magnitudes (17.5<I<24), at any given redshift we can 
only observe galaxies in a specific, redshift-dependent, ab- 
solute magnitude range. It is usual to describe the sample 
radial incompleteness by defining the selection function in 
terms of the galaxy luminosity function <p(A4) 



f M ^<p(M)dM ' 



S(r,M c ) = 



(10) 



Here Mbir) and Mf{r) are the B-band absolute mag- 
nitudes which correspond, at distance r, to the I-band lim- 
iting apparent magnitudes mj, = 17.5 and nif — 24 respec- 
tively (see discussion at the end of the section). Since eq. 

also depends on galaxy colors, we compute its mean 
value at distance r by a weighted average over the popu- 
lation mix observed at distance r. 

The VVDS luminosity function (LF) has been derived 
in Paper II and is characterized by a substantial degree of 
evolution over the redshift range 0< z<1.5. Therefore, we 
estimate tp(M.) in the B band at any given position in the 
redshift interval [0, 1.5] by interpolating, with a low order 
polynomial function, the Schechter shape parameters a 
and Ai* given in Table 1 of Paper II. 

Assuming M. C B = —15 in eq. which corresponds to 
the limiting absolute magnitude over which the LF of the 
VVDS-02h sample can be robustly constrained in the low- 
est redshift bins, the selection function exponentially falls 
by nearly 2 orders of magnitude in the redshift range up 
to z<1.5. Thus, the density field reconstruction strongly 
depends on the radial selection function used especially 
at high redshifts, where eq. EI] can be affected by possi- 
ble systematics in the determination of the LF or in the 
measurements of faint magnitudes. Therefore, we will also 
analyze volume-limited sub-samples, which are essentially 
free from these systematics. 

Since in a magnitude limited survey progressively 
brighter galaxies are selected as a function of redshift, a 
volume-limited sample also allows us to disentangle spu- 
rious luminosity-dependent effects from the measurement 
of the redshift evolution of the biasing function. 

3.3. The Redshift Sampling Rate 

As for most redshift surveys, the VVDS does not target 
spectroscopically all the galaxies which satisfy the given 
flux limit criteria in the selected field of view (see Paper I) . 
Because of the sparse sampling strategy, we have to correct 
the density estimator with a sampling rate weight $ 2 in 
order to reconstruct the real underlying galaxy density 
field in a statistically unbiased way. 

The VVDS redshift sampling rate is the combination 
of two effects: i) only a fraction of the galaxies (~ 40%, 
see Paper I) in the photometric sample is targeted (target 
sampling rate), ii) and only a fraction of the targeted ob- 
jects (~ 80% see Paper I) yield a redshift (spectroscopic 
success rate). We can model this correction term, by as- 
suming, in first approximation, that the sampling rate de- 
pends only on the apparent magnitude. Since the VVDS 
targeting strategy is optimized to maximize the number of 
slits on the sky, the selection of faint objects is systemati- 
cally favored. Inversely, the ability of measuring a redshift 
degrades progressively towards fainter magnitudes, i.e., for 
spectra having lower signal-to-noise ratios (the spectro- 
scopic success rate decreases from > 90% at 1^22 down 
to ~ 60% at I~24.) These two opposite effects conspire to 
produce the magnitude-dependent sampling rate function 
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Fig. 1. The VVDS redshift sampling rate in the four- 
passes VVDS-02h-4 region is plotted versus the observed 
apparent magnitude in the I-band. The mean redshift 
sampling rate is ~ 0.3. 

shown in Fig. ^ Clearly, with such an approximation, we 
neglect any possible dependence of the sampling rate from 
other important parameters such as for example surface 
brightness, spectral type or redshift. However, in in §3.2 of 
Paper II we showed, using photometric redshifts, that any 
systematic sampling bias introduced by a possible redshift 
dependence of the spectroscopic success rate is expected 
to affect only the tails of our observed redshift distribution 
(z<0.5 and z>1.5) i.e., redshift intervals not considered 
in this study (see §3). 

We describe the VVDS sampling completeness $ z (rn) 
at a given magnitude m, as the fraction of objects with 
measured redshifts N z over the number N of objects de- 
tected in the photometric catalog. 



J2i=x w ( m - m »> dm) 



$*(m) 



with the window function w defined as 
w(m — mi, dm) - 



1 if \m — mi\ < dm/2 
otherwise 



(11) 



(12) 



Working in low resolution mode (i.e allowing up to 4 
galaxies to have spectra aligned along the same disper- 
sion direction with a typical sky separation of 2 arcmin) 
and observing the same region of sky 4 times for a to- 
tal of about 16 hours (four-passes strategy), we can mea- 
sure redshifts for nearly 30% of the parent photometric 
population, as shown in Fig. ^ in other terms, on aver- 
age, nearly one over three galaxies with magnitude I<24 



has a measured redshift in the four-passes VVDS region. 
This high spatial sampling rate is a critical factor for 
minimizing biases in the reconstruction of the 3D den- 
sity field of galaxies. In particular we note that our I<24, 
z<1.5 VVDS-02h-4 sample is characterized by an effec- 
tive mean inter-particle separation in the redshift range 
[0,1.5] ((r) ~ 5.1 /i _1 Mpc ) which is smaller than that 
of the original CFA sa mple ((r) ~ 5.5/i _1 Mpc) used by 
iDavis fc Huchral 1 198ll) to reconstruct the 3D density field 
of the local Universe (i.e., out to ~ 80 /i _1 Mpc ). At 
the median depth of the survey, i.e., in the redshift in- 
terval 0.75±0.05, the mean inter-particle separation is 4.4 
/i _1 Mpc, a value nearly equal to the mean inter-particle 
separation at the median depth of the 2dFGRS. We fi- 
nally note that in the redshift range [0.7,1.35], which is 
also covered by the DEEP2 survey, the VVDS mean-inter- 
particle separation is 5.5 /i _1 Mpc compared to the value 
~ 6.5 /i _1 Mpc inferred from the values quoted by Coil et 
al. (2004) for their most complete field, currently covering 
0.32 deg 2 . 

By dividing the VVDS-02h-4 field in smaller cells and 
repeating the analysis, we conclude that the sampling rate 
does not show appreciable variations, i.e., the angular se- 
lection function can be considered constant for the pur- 
poses of our analysis. This corresponds to the fact that 
the success rate in redshift measurement in each VIMOS 
quadrant (i.e., the spectroscopic success rate per mask) is, 
in good approximation, constant and equal to ~80% (see 
Paper I). 



3.4. Shot Noise 

In a flux-limited sample, the shot noise in the density field 
is an increasing function of distance (see eq. (jBJl). One 
may correct for the increase of the mean VVDS inter- 
particle separation as a function of redshift (and thus 
the increase of the variance of the density field) by op- 
portu nely increasing the length of the smoothing window 
(e.g., IStrauss fc Willickl Il995l . However, since we are in- 
terested in comparing the fluctuations recovered on the 
same scale at different redshifts in a flux limited survey, 
we take into account the decreased sampling sensitivity of 
the survey at high redshift in an alternative way. 

We deconvolve the signature of this noise from the 
density maps by a pplying the Wiener filter technique (cf. 
IPress et all Il992h which provides the minimum variance 
reconstruction of the smoothed density field, given the 
map of the noise and the a pri ori knowledge of th e un- 
derlying power spectrum (e.g., lLahav et all 1199^1 . The 
application of the Wiener denoising procedure to the spe- 
cific geometry of the VVDS sample is described in detail 
in appendix A. 

Here we note that the Wiener filter requires a model 
for the underlying 3D power spectrum P(k,z) which we 
compute, over the frequency scales where the correlation 
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Fig. 2. The real- and redshift-space rms fluctuations of the 
flux-limited VVDS sample recovered using eq. and 
the results of the correlation function analysis presented 
in Paper III are plotted at six different redshifts in the 
interval 0.4< z<1.7. 

function of VVDS galaxies is well constrained (0.06 < k < 
10), as (see eq. in Appendix A): 



P(k,z) = 4tt 



ro{z 



fc3- 7 (z) 



■T(2 - 7(2)) sin 



(13) 



where the normalization ro(z) and the slope 7(2) of 
the correlation function at redshift z have been derived 
by interpolating the values measured in various reds hift 
interva ls of the VVDS-02h volume by iLe Fevre et alJ 
(2005b) (hereafter Paper III). 

The variance of the galaxy distribution on a 8 h~ 1 Mpc 
scale in the VVDS-02h sample can be obtained by inte- 
grating eq. (|13fl using a TH window of radius 8 h~ 1 Mpc 
and the (^0,7) parameters of the VVDS correlation func- 
tion 



1 

2^2 



P{k)F£k 2 dk 



(14) 



where Fk is the Fourier transform of the TH filter (see eq. 
I|45|) in Appendix A). Note that, by integrating the power 
spectrum down to k — * 0, i.e., extrapolating the power 
law shape of eq. (|T3J) beyond L ~ 100 /i^Mpc (k ^0.06), 
one would revise upwards the value of as by ~ 2% (at 
z=0.35) and by ~ 4% at z—lA. Since, however, the am- 
plitude of the power spectrum on large scales is expected 
to downturn and to be systematically lower than predicted 
by eq. I|13|) . we safely conclude that, with our computation 
scheme, the inferred cr 8 value should be biased low by no 



more than ~ 2% and 4% in the first and last redshift bin, 
respectively. 

Projecting the results from real-space into the redshift- 
distorted space (see §5), i.e., implementing the effects 
of large-scale streaming motions, we infer that the rms 
of galaxy fluctuations are <7s = [0.55±0.10, 0.62±0.12, 
0.71±0.11, 0.69±0.14, 0.75±0.14, 0.92±0.20] at redshift 
z=[0. 35,0.6,0.8, 1.2, 1.2, 1.65] (see Fig. 0). Thus, the am- 
plitude of as for a flux- limited I<24 sample increases as 
a function of redshift by nearly 70% between z~0.3 and 
z~1.7. 

3.5. The 3-Dimensional VVDS Density Field 

The VVDS-02h galaxy density field reconstructed on a 
scale R= 5 h~ 1 Mpc and in the redshift bin 0.8< z<l.l is 
visually displayed in the left most panel of Fig. Note 
that the chosen smoothing scale nearly corresponds to the 
mean inter-galaxy separation in the VVDS-02h sample at 
1. Fig. shows the regular patterns traced by over- 
and under-dense regions in the selected redshift interval. 
Specifically, we note that, in this redshift slice, there are 
over- and under-dense regions which extend over charac- 
teristic scales as large as ~ 100 ft, _1 Mpc. A more complete 
discussion of the "cartog raphy" in such deep reg ions of the 
Universe is presented bv lLe Fevre et all l 2005eh . 

In the same figure, we also display the density 
field reconstructed in an analogous redshift range, us- 
ing the GA LICS simulat i on (G alaxies in Cosmological 
Simulations. lHatton et all l|2003j) '). GALICS combines cos- 
mological simulations of dark matter with semi-analytic 
prescriptions for galaxy formation to produce a fully re- 
alistic deep galaxy sample. In particular we plot the den- 
sity field of the I<24 flux- limited simulation as well as 
the density field recovered after applying to the pure flux- 
limited simulation all the VVDS target selection criteria 
(see §4.1). No qualitative difference between the density 
fields reconstructed before and after applying to the sim- 
ulation all the survey systematics is seen. 

Clearly, a more quantitative assessment of the robust- 
ness and reliability of the VVDS overdensity field can be 
done by studying its PDF. 

4. The PDF of Galaxy Fluctuations 

Once the three-dimensional field of galaxy density con- 
trasts S g has been reconstructed on a given scale R, one 
can fully describe its properties by deriving the associated 
PDF g{8g). This statistical quantity represents the nor- 
malized probability of having a density fluctuation in the 
range (5 g , S g + dS g ) within a region of characteristic length 
R randomly located in the survey volume. 

While the shape of the PDF of mass fluctuations at 
any given cosmic epoch is theoretically well constrained 
from CDM simulations (see next section), little is known 
about the observational PDF of the general population of 
galaxies in the high redshift Universe. Even locally, this 
fundamental statistics has been often overlooked (but see 
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Fig. 3. 3D density field traced by the galaxy distribution in the VVDS-02h Field (left, 1641 galaxies), in the flux- 
limited (I<24) GALICS simulation (right, 9450 galaxies), and in the flux-limited GALICS sample after applying the 
VVDS target selection criteria (center, 1656 galaxies). Data span the redshift interval [0.8,1.1]. In each cone, the galaxy 
distribution is continuously smoothed using a TH window function with R=5 ft, _1 Mpc which nearly corresponds to 
the mean inter-particle separation in this same redshift interval. The cone metric was computed assuming a ACDM 
cosmology and the correct axis ratio between transversal and radial dimensions has been preserved. The cones have 
approximate transverse dimensions of 28 h~ 1 Mpc at z—1 and extend over 527 ft, Mpc. 



lOstriker et alJ |2003)). Notwithstanding, the shape of the 
galaxy PDF is strongly sensitive to the effects of gravita- 
tional instability and galaxy biasing, and its redshift de- 
pendence encodes valuable information about the origin 
and evolution of galaxy density fluctuations. 

The shape of the PDF can be characterized in terms 
of its statistical moments. In particular the variance of 
a zero-mean field (such as the overdensity field we are 
considering) is 



4 = (S 2 g )R = / S 2 g g R (6 g )dS g . (15) 
Jo 

Hig her moments can be straightforwardly derived (see 
iBernardeau et aill2002l for a review). In the following, we 
will drop the suffix R, unless we need to emphasize it. 



4.1. Estimating Reconstruction Systematics Using 
Mock Catalogs 

For the purposes of our analysis, it is imperative to check 
that the various instrumental selection effects as well as 
the VVDS observing strategy are not compromising the 
determination of the PDF of galaxy fluctuations. In this 
section, we explore the region of the parameter space (es- 
sentially redshift and smoothing scales) where the PDF of 
VVDS overdensities traces in a statistically unbiased way 
the underlying parent distribution. 

Possible systematics can be hidden in the recon- 
structed PDF essentially because i) the VVDS redshift 
sampling rate is not unity, ii) the slitlets are allocated 
on the VIMOS masks with different constraints along the 
dispersion and the spatial axis, and iii) the VIMOS field 
of view is splitted in four different rectangular quadrants 
separated by a vignetting cross. 
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We have addressed point iii) by designing a specific 
telescope pointing strategy which allows us to cover in 
a uniform way the survey sky region (see the telescope 
pointing strategy shown in Fig 1 of Paper I). With the 
adopted survey strategy, we give to each galaxy in the 
VVDS-02h-4 field four chances to be targeted by VIMOS, 
thus increasing the survey sampling rate (nearly 1 galaxy 
over 3 with magnitude I<24 has a measured redshift.) 

Concerns about points i) and ii) can be directly ad- 
dressed using galaxy simulations covering a cosmolog- 
ical volume comparable to the VVDS one. Thanks to 
the implementation of the Mock Map Facility (MoMaF 
iBlaizot et all 12003]) . it is possible to convert the GALICS 
3D mocks catalog into 2D sky images, and handle the 2D 
projection of the simulation as a pseudo-real imaging sur- 
vey. Polio et al. (2005), have then built a set of 50 fully 
realistic mock VVDs surveys from the GALICS simula- 
tions to which the whole observational pipeline and bi- 
ases has been applied. By comparing specific properties of 
the resulting pseudo-VVDS sample with the true under- 
lying properties of the pseudo-real Universe from which 
the sample is extracted, we can directly explore the ro- 
bustness, as well as the limits, of the particular statistical 
quantities we are interested in. 

In brief these include addeding to the 3D galaxy mocks 
a randomly simulated distribution of stars to mimic the 
same star contamination affecting our survey. Next, we 
have masked the sky mocks using the VVDS photomet- 
ric masks, i.e., we have implemented the same geomet- 
rical pattern of excluded regions with which we avoid 
to survey sky regions contaminated by the presence of 
bright stars or photometric defects. Then, we have ex- 
tracted the spectroscopi c targets by applyin g the target 
selection code (VMMPS, iBottini et all l2005|) to the sim- 
ulated 2D sky distribution. To each GALICS redshift, 
which incorporates the Doppler contribution due to galaxy 
peculiar velocities, we have added a random component 
to take into account errors in z measurements. Finally, 
we have processed the selected objects implementing the 
same magnitude-distribution of failures in redshift mea- 
surements which characterizes the first-epoch observations 
of the VVDS survey (see Paper I and FigQJ. 

Since GALICS galaxies have magnitudes simulated in 
the same 4 bands surveyed by VVDS (B,V,R,I), we have 
applied the K-correction to obtain rest frame absolute 
magnitudes and we have empirically re-derived all the se- 
lection functions for the mock catalogs according to the 
scheme presented in §3. In this way we can also check the 
robustness of the techniques we apply for computing ab- 
solute magnitudes and selection functions (see Paper II). 

The PDF of galaxy overdensities obtained from the 
mock samples has been finally compared to the PDF of 
the parent population. For brevity, in what follows, we 
will call s-samples (survey-samples) the mocks simulating 
the VVDS redshift survey and p-sample (parent sample) 
the whole GALICS simulation flux-limited at 17.5<I<24. 

The density contrasts have been calculated as de- 
scribed in §3. In the following, we will restrict our analysis 
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Fig. 4. Tests of the PDF reconstruction scheme using 
the mocks VVDS samples extracted from GALICS. The 
differential (G(y g )) probability distribution functions of 
y g = l + Sg for the "observed" s-sample (dotted line, shad- 
owed histogram) , and for the parent p-sample (solid- line) . 
Note that the plotted histograms actually corresponds 
to G{y g ) = hi(l0)y g g(y g )) since the binning is done in 
log(y g ). The logarithmic PDFs are computed for density 
fields smoothed using TH filters of different sizes (indi- 
cated on the top of each panel). 



to the set of smoothing scales in the interval R =(5,10) 
/i _1 Mpe. The choice of these particular limits is motivated 
by the fact that 5 h~ 1 Mpc is the minimum smoothing scale 
for which the reconstructed density field is unbiased over a 
substantial redshift interval (see discussion below). Note, 
also, that below this typical scale, linear regimes approxi- 
mations, largely used in this paper, do not hold anymore. 
The upper boundary is constrained by the transverse co- 
moving dimensions covered by the first-epoch VVDS data 
(see §2), which is still too small for being partitioned using 
bigger scale-lengths without introducing significant noise 
in the reconstructed PDF (see the transverse comoving 
dimension of the VVDS-02h field quoted in §2). 
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Fig. 5. The cumulative distribution function of density 
contrasts y g = 1 + 5 g , on scales R=5 /i _1 Mpc, as recov- 
ered in different redshift intervals for the s-sample (dotted 
line), and for the parent p-sample (solid- line). The flat- 
tened pedestal at the low-density end of the cumulative 
distribution is due to low density regions in the p-sample 
that are spuriously sampled as empty regions (S g = — 1) 
in the s-sample. 



low density regions of the Universe (log(l + S g ) <; —0.5) 
when the galaxy distribution is smoothed on scales of 5 
/i _1 Mpc. Fig.Elshows that incompleteness in underdense 
regions is a function of redshift, with the bias in the low- 
density tail of the PDF developing and increasing as the 
redshift increases. As a rule of thumb, the PDF of the s- 
sample starts to deviate significantly from the parent PDF 
when the mean inter-galactic separation (r) of the survey 
sample is larger than the scale R on which the field is 
reconstructed. Since we measure (r(z — 1)) ~ 5 /i _1 Mpc 
((r(z = 1.5)) - 8 h^Mpc ), the PDF of the density 
field recovered using a TH filter of radius 5 ft^Mpc is 
effectively unbiased (at least over the density range we 
are interested in, i.e., log{\ + 6) > —1) only if the sample 
is limited at z<; 1. Therefore, in the following, results 
obtained for R=5 h~ 1 Mpc are quoted only up to z=l. 

On scales R> 8/i _1 Mpc , the agreement between the 
PDFs of s- and p-samples holds true also for volume- 
limited subsamples Specifically, the 2nd and 3rd mo- 
ments of the PDF of overdensities recovered using volume- 
limited s-samples on these scales are within la of the 
corresponding values computed for the parent, volume- 
limited, p-samples in each redshift interval of interest up 
to z = 1.5. 

To summarize, the results of simulated VVDS obser- 
vations presented in this section show that, at least on 
scales R>8 /i^Mpc, the VVDS PDF describes in an un- 
biased way the general distribution properties of a sample 
of 1=24 flux- limited galaxies up to redshift 1.5. In other 
terms, the VVDS density field sampled in this way is es- 
sentially free from selection systematics in both low- and 
high-density regions, and can be meaningfully used to infer 
the physical bias in the distribution between galaxy and 
matter. Obviously, the representativeness of the measured 
PDF of VVDS overdensities with respect to the "univer- 
sal" PDF up to 2 = 1.5 is a different question. Since the 
volume probed is still restricted to a limited region of space 
in one field, the shape and moments of the galaxy PDF 
derived from the VVDS-02h are expected to deviate from 
the "universal" PDF of galaxies at this redshift because 
of cosmic variance. Reducing the cosmic variance is one of 
the main goals of extending the VVDS to 4 independent 
fields. Anyway, our 50 mock realisations of the VVDS-02h 
sample allow us to estimate realistic errors that include 
the contribution from cosmic variance. 



The PDF of the galaxy density contrasts computed 
using the s-sample is compared to the parent distribution 
inferred from the p-sample in Fig. We conclude that 
the distribution of galaxy overdensities of the s-samples 
for R = 8 and 10 ft^Mpc scales is not biased with re- 
spect to the underlying distribution of p-sample galaxy 
fluctuations. Thus, the VVDS density field reconstructed 
on these scales is not affected by the specific VVDS ob- 
servational strategy. 

It is evident in Fig.^Jthat the VVDS redshift sampling 
rate is not high enough to map in an unbiased way the 



4.2. The PDF of VVDS Galaxies: Results 

Let us then investigate the evolution, as a function of the 
lookback time, of the observed PDF of VVDS galaxy fluc- 
tuations. 

In an apparent magnitude-limited survey such as the 
VVDS, only brighter galaxies populate the most distant 
redshift bins, whereas fainter galaxies are visible only at 
lower redshifts. As more luminous galax ies tend to clus- 
ter more strongly th an fainter ones (e.g., lHamiltorl [l988: 
ICroton et alll2004L the PDF of galaxy fluctuations is ex- 
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Fig. 6. Top: The PDF per units of galaxy overdensities (G(y) = ln(10)yg(yj) is plotted for the volume limited VVDS 
sample (M% = —20 + 5 log ft,) at three different rcdshifts. The PDFs are computed for density fields recovered using 
TH filters of size 8 and 10 fr^Mpc. Bottom: the corresponding cumulative distributions. Errorbars (which, for clarity, 
are plotted only for the redshift bin 1.25 < z < 1.5) represent the Poissonian uncertainties. The flattened pedestal at 
the low-density end of the cumulative distribution is contributed by empty regions (S g = —1). 



pected to be systematically biased as a function of red- 
shift. 

This effect is clearly seen in the first correlation anal- 
ysis of the VVDS (paper III), and can be minimized by 
selecting a volume-limited sample. Therefore, we have de- 
fined a subsample with absolute magnitude brighter than 
M% = -20 + 51ogfr in the rest frame B band (~ 1350 
galaxies with 0.7< z<l.5 in the VVDS-02h field, - 800 of 
which are in the VVDS-02h-4 region). 



This threshold corresponds to the faintest galaxy lu- 
minosity visible at redshift z=1.5 in a 1=24 flux-limited 
survey and it is roughly 0.6 magnitudes brightcr(fainter) 
than the value of Ai* B recovered at z~0(~1.5) using the 
VVDS data (see Paper II). The median absolute magni- 
tude for this volume-limited sample is ~ —20.4. 

We note, however, that the populations of galaxies 
with the same luminosity at different redshifts may actu- 
ally be different. As we have shown in Paper II, we mea- 
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sure a substantial degree of evolution in the luminosity of 
galaxies, and, as a consequence, with our absolute mag- 
nitude cut we are selecting M* + 0.6 galaxies at 2=1.5, 
but A4* — 0.6 galaxies at z=0. Thus, the clustering signal 
at progressively earlier epochs may not be contributed by 
the progenitors of the galaxies that are sampled at later 
times in the same luminosity interval. 

The PDF of density fluctuations, in various redshift 
intervals, and traced on scales of 8 and 10 ft, _1 Mpc by 
VVDS galaxies brighter than A4%, is presented in Fig. EI 
Note that the analysis of the previous section guarantees 
that, on these scales, the VVDS PDF fairly represents the 
PDF of the real underlying population of galaxies up to 
.z=1.5. Fig.EJshows how the shape of the measured galaxy 
PDF changes across different cosmic epochs. 

A Kolmogorov-Smirnovtest confirms that the PDFs at 
different cosmic epochs are statistically different (i.e., the 
null hypothesis that the three distributions are drawn from 
the same parent distribution is rejected with a confidence 
P fcs > 1-10- 6 ). 

In particular the peak of the PDF in the lowest redshift 
interval is shifted towards smaller values of the density 
contrast 6 g when compared to the peak of the PDF in the 
highest redshift bin. Moreover, the shape of the distribu- 
tion, also shows a systematic "deformation" . Specifically, 
we observe the development of a low-5 tail in the PDF as a 
function of time on both scales investigated. In other terms 
the probability of having low density regions increases as a 
function of cosmic time. For example, underdense regions, 
defined as the regions where the galaxy density field is 
log(l + S g ) < —0.5 on a R= 8 /i _1 Mpc scale, occupy a 
fraction of nearly 35% of the VVDS volume at redshift 
0.7< 2<1.0, but only a fraction of about 25% at redshift 
1.25<z<1.5. Similar trends are observed when lowering 
the absolute luminosity threshold of the volume limited 
sample (and consequently lowering the upper limit of the 
redshift interval probed) or when modifying the binning 
in redshift space. 

If galaxies are faithful and unbiased tracers of the un- 
derlying dark matter field, then both this effects, the peak 
shift and the development of a low density tail may be 
qualitatively interpreted as a direct supporting evidence 
for the paradigm of the evolution of gravitational cluster- 
ing in an expanding Universe. At variance with overdense 
regions which collapse, a net density deficit (5 < 0) in 
an expanding Universe brings about a sign reversal of the 
effective gravitational force: a density depression is a re- 
gion that induces an effective repulsive peculiar gravity 
(Peebles 1980). If gravity is the engine which drives clus- 
tering in an expanding Universe, we thus expect that, as 
time goes by, low density regions propagate outwards and 
a progressively higher portion of the cosmological volume 
becomes underdense. 

The observed evolution in the PDF could also indicate 
the presence of a timc-dcpcndcnt biasing between matter 
and galaxies. As a matter of fact, it can be easily shown 
that a monotonic bias, increasing with redshift, offers a 
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Fig. 7. Redshift evolution of the standard deviation (up- 
per panel and of the skewness (lower panel) of the galaxy 
PDF on a scale R=8 ft, _1 Mpc for galaxies brighter than 
M = —20 + 51og/i. The corresponding local values, es- 
timated on the same scale by Croton et al. 2004 using a 
subsample of the 2dFGRS having nearly the same median 
absolute luminosity of our sample, are represented with 
triangles. Error bars represent la errors, and, in the case of 
VVDS measurements, include the contribution from cos- 
mic variance. The errorbar on as of the 2dFGRS is smaller 
than the symbol size. 

natural mechanism to re-map the galaxy PDF into pro- 
gressively higher intervals of density contrasts. 

We can better discriminate the physical origin of the 
observed trends, i.e., if they are purely induced by grav- 
itation or strengthened by the collateral and cooperative 
action of biasing, by studying the evolution of the PDF 
moments. In Fig. the redshift evolution of the rms (a) 
and skewness (S3) of the overdensity fields (see table 1) for 
the M.b < —20 + 5 log h sample are shown and compared 
to local measurements. 

Following the standard convention within the hier- 
archical clustering model, we define the skewness S3 in 
terms of the volume-averaged two- and three-point cor- 
relation functions (S3 = £,3/ £,2) noting that in the case 
of a continuous 5 g -field with zero mean this expression 
reduces to S3 = (Sg)/(Sg) 2 . We do not derive the mo- 
ments < S 2 > and < 6g > of the PDF by directly ap- 
plying the computation scheme given in eq. H15[). but by 
correcting the count-in-cells statistics for discret eness ef- 
fects using the P oissonian shot-noise model (e.g., IPeeblesl 
Il98fi lFrvLfl985L cfr. eqs. (374) and (375) of Bernardeau 
et al. 2002, possible biases introduced by this estimation 
technique are discussed by Hui & Gaztafiaga 1999.). The 
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corresponding values of a and S3 for the local Universe (in 
redshift space) have been derived by Croton et al. 2004 us- 
ing the 2dFGRS. Here we plot the values corresponding 
to their —21 < M — b \ogh < — 20 subsample, which actu- 
ally brackets the median luminosity of our volume limited 
sample. 

We can see that the rms amplitude of fluctuations 
of the VVDS density field, on scales 8 /i~ 1 Mpc, is with 
good approximation constant over the full redshift base- 
line investigated, with a mean value of 0.94±0.07 over 
0.7<z<1.5. 

While the strength of clustering of galaxies brighter 
than M. < —20 + 51og/i does not change much in this 
redshift interval, each VVDS measurement is lower than 
the value inferred at z~0 by Croton et al. 2004. In partic- 
ular, our mean value is ~10% smaller than the 2dFGRS 
value and the difference is significant at ~ 2a level. 

The skewness S3 , which measures the tendency of grav- 
itational clustering to create asymmetries between under- 
dense and overdense regions, decreases as a function of 
redshift. We observe a systematic decrement not only in- 
ternally to the VVDS sample, but also when we compare 
our measurements with the z—0 estimate. This trend is 
caused by the development of the low-<5 tail in the PDF as 
a function of time on both the R=8,10 h~ 1 Mpc scales and 
reflects the fact that the probability of having underdense 
regions is greater at present epoch than it was at z~1.5 
(where its measured value is ~ 2a lower.) 

The amplitudes of the rms and skewness of galaxy 
overdensities show an evolutionary trend dissimilar from 
that predicted in first and second order perturbation the- 
ory for the gravitational growth of dark matter fluc- 
tuations (see Bernardeau et al. 2002 for a review). 
According to linear perturbation theory the amplitude 
of the rms of mass fluctuations scales with redshift as 
in eq. (|18fl while second order perturbation theory pre- 
dicts that, on the scales where the quasi-linear approx- 
imation holds, the growth rate of < Sg > and vari- 
ance < 6 2 > 2 are syncronized so that the skewness S3 
of an in i tially Gaussian fields should remain constant 
llPeeblesl Il980t lJuszkiewicz. Bouchet fc Colombit Il993l 
iBernardeaul Il993(l . 1 Furthermore, in Le Fevre et al. 
(2005c) we show that even the general shape of the galaxy 
PDF deviates from a lognormal distribution, i.e., from the 
profile in terms of which the mass PDF is generally ap- 
proximated (see §5). Therefore, we conclude that the PDF 
evolution is not caused by gravity alone; the redshift scal- 
ing of its global shape and moments effectively indicates 
the presence of a time evolving bias. 

We can deconvolve the purely gravitational signature 
and investigate properties and characteristics of the bias- 
ing between matter and galaxies by comparing the galaxy 
PDF to the corresponding statistics computed for mass 



1 Note that the observed redshift evolution of the skewness 
is just the opposite of what is expected also in generic dimen- 
sional non-Gaussian models where 5*3 is predicted to increase 
with redshift. 



fluctuations. Thus, we now turn to the problem of deriv- 
ing the PDF of mass fluctuations. 

5. The PDF of Mass Fluctuations in 

Redshift- Distorted Comoving Coordinates 

The VVDS survey is providing a rich body of redshift data 
for mapping the galaxy density field in extended regions 
of space and over a wide interval of cosmic epochs. On the 
contrary, the direct determination of the underlying mass 
density field and its associated PDF is a less straightfor- 
ward process. Nonetheless we may gain insight into the 
mass statistics by using simulations and theoretical argu- 
ments. 

In the standard picture of gravitational instability, the 
PDF of the primordial cosmological mass density fluctu- 
ations is assumed to obey a random Gaussian distribu- 
tion. Once the density fluctuations reach the non-linear 
stage, their PDF significantly deviates from the initial 
Gaussian profile and a variety of phenomenologic al mod- 
els have been proposed to describe its shape (e.g., ISaslawL 
1198,4 lLahav et all Il99.^ . In particular, it is well estab- 
lished in CDM models that when structure formation has 
reached the nonlinear regime, the density contrasts in co- 
moving space f(S) fo llow, to a good approximation, a log- 
normal distribution llColes fc JonesL Il99lt iKofman et a 1 



Il99l iTavlor fc WattsL 1200(1 iKavo. Tarn va fc Sutol. bOOlli. 



IW = mr ex P 



{ 



[ln(l + 6) + m 2 /2} 2 
2uj 2 



} (16) 



This approximation becomes poor in th e highly 
non-linear regim e feernardeau fc KofmanL 1199.4 
lUeda fc Yokovamal Il996t e.gl. The PDF of mass 
overdensities f(S) is characterized by a single parameter 
(w) that is related to the variance of the 5-field as 



lu 2 = ln[l + (S 2 )} 



(17) 



At high redshifts, the variance a^ over sufficiently 
large scales R (those explored in this paper) may be easily 
derived using the linear theory approximation: 



<tr{z) = <tr(z = 0)D{z) 



(18) 



where D(z) is the linear grow th rate of density fluctuation s 
normalized to unity at z=0 l|Heatl Il977t iHamiltonl . l200lj) . 

The lognormal approximation formally describes the 
distribution of matter fluctuations computed in real co- 
moving coordinates. On the contrary, the PDF of galax- 
ies is observationally derived in redshift space. In or- 
der to map properly the mass overdensities into galaxy 
overdensities the mass and galaxy PDFs must be com- 
pu ted in a common reference fram e. It has been shown 
bv lSigad. Branchini fc Dekell l)2000|) that an optimal strat- 
egy to derive galaxy biasing is to compare both mass and 
galaxy density fields directly in redshift space. Implicit in 
this approach is the assumption that mass and galaxies 
are statistically affected in the same way by gravitational 
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perturbations, and thus, that there is no velocity bias in 
the motion of the two components. 

A general model which allows the explicit computa- 
tion of the statistical distortions caused by peculiar ve- 
locities has been proposed by Kaiser (1987). This applies 
in the linear regime (i.e., on large scales) and in the lo- 
cal Universe where redshift and distances are linearly re- 
lated. At cosmological distances z, however, the mapping 
between real comoving coordinates (x) and redshift co- 
moving coordinates (y), i.e., the pseudo-comoving coordi- 
nates inferred on the basis of the observed redshifts, is less 
trivial, and we proceed to obtain it in the following. 

In an inhomogeneous Universe, galaxies have 
motions above and b e yond t heir Hubble v elocity 
(e.g.. iGiovanelli et all Il998t iMarinoni etafl Il998t 
iBranchini et all l200l[) . As a consequence, Doppler spec- 
tral shifts add to the cosmological signal and the observed 
redshift (if) is given by 



z = z + 



m (1 




(19) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



where z is the cosmological redshift in a uniform 
Friedman-Robertson- Walker metric and where J7(x) = 
v(x) ■ f = |v(x)|/x is the radial component of the peculiar 
velocity {\x is the cosine of the angle between the peculiar 
velocity vector and the line-of-sight versor f). 

The redshift comoving distance of a galaxy at the ob- 
served redshift z is thus 



y 

where 



c 

Hq jo 



2+2(1 + 2) 



1 



E(x) 



(20) 



E(z) = [fi m (l + zf + (l-n m -fi A )(l + zf +ft A ] 1/2 .(21) 
At high redshifts (z > — ), we can write 



V = 



J—dx + -{l + z)E{z)- 1 
o E {x) c 



(22) 



which, in turns, gives the coordinate transformation from 
real comoving space x to the redshift comoving space y 



y 



l+p{z) 



<7(x) 



In this mapping, the cosmological term p{z) , 



P(*) = 



1 + z 
H E(z) 



(23) 



(24) 



is a correcting factor which takes into account the fact 
that, at high redshifts, distances do not scale linearly 
with redshift, and, thus, that peculiar velocities cannot 
be simply added to redshift space positions as in the local 
Universe. 

The galaxy density field in the redshift-distorted space 
is related to the galaxy density in real space by the 



Fig. 8. Redshift scaling of the rms mass fluctuations in 
sphere of 8 h~ 1 Mpc radius. Diamonds represent as as com- 
puted from the ACDM Hubble volume simulation in real 
comoving space (x-space), while triangles represent the 
corresponding values recovered in the redshift perturbed 
comoving coordinates (y-space). The solid line is the an- 
alytical prediction for the scaling of as in the y-space ob- 
tained using eq. (|30[1 . while the dotted line represents the 
x-space evolution predicted in real space. In both cases the 
power spectrum of perturbations is the same and has been 
normalized in order to match the simulation specifications 
(a x (z = 0) = 0.9). 

Jacobian of the transformation between the two coordi- 
nate systems 



Pv(y) = Px(x) l+p(z) 



l+p(z) 



dU(x) 



dx 



(25) 



At sufficiently large distances from the observer, ne- 
glecting the survey selection function (i.e., considering a 
volume-limited redshift survey) and at first order in per- 
turbations we obtain 



<5y(y) = ^( x ) -p( z ) 



dU(x) 



(26) 



The second term on the right hand side can be evalu- 
ated using linear-regime approximations and gravitational 
instability theory. In comoving coordinates it is given by 



dU(x) 



M 2 f(z)H(z) 
1 + z 



(27) 



where f = dhiD/dlna is the logarithmic derivative of the 
linear growth rate of density fluctuations with respect to 
the expansion factor a(t). At redshift z (corresponding to 
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Fig. 9. One-point PDFs of dark matter fluctuations 
(shaded area) computed using the Hubble volume ACDM 
cosmological simulation in 4 different redshift ranges over 
a volume which mimics the geometry of the VVDS sam- 
ple. The mass PDF has been recovered in the redshift 
comoving space by smoothing the mass-particle distribu- 
tion with a TH window of size R=8 h~ 1 Mpc. Note that 
the plotted histogram actually corresponds to F(y) = 
(lnlO)yf(y) because the binning (dlogy=0.1) is done in 
log(y) = log(l + S). The dotted line represents the log- 
normal approximation derived in the real comoving space 
using eq. (|16() . The solid curve represents the lognormal 
approximation computed adopting the variance parameter 
lo (shown in the inset) theoretically inferred using equa- 
tion l|30[) . which models peculiar velocity distortions as a 
function of redshift. 

the comoving position x) a useful approximation is given 
by: 



f(z) ~ il 3 n ( 5 E(zy 6 ^{l + z) 9 / 5 (28) 

fsee lMartell dl99lh : lLahav et alJ jl99lln . 

By combining the previous results we obtain 



S y (y) = 4(x)[l + /i 2 f(z)]. 



(29) 



which reduces to the iKaiserl Jl987l) correction when 2=0 
The relation between the azimuthally averaged vari- 
ances measured in real and redshift comoving space is 



cjV{z)= l + lf( z ) + ^f 2 (z) V V(z) 



(30) 



We have tested the validity of eq. I|30|) in the high 
redshift domain using the Hubble volume N-body simula- 
tions carried out by the Virgo consortium ( Colberg et all 



l2000h . This is a large numerical experiment which allows 
the simulation of mass surveys along the observer past 
light cone. The simulated mass distribution is computed in 
a ACDM cosmogony with parameters f2 m = 0.3, £1a = 0.7 
and Hq — 70 km s _1 Mpc _1 . The volume covered by this 
N-body simulation is large enough that the mass survey 
extracted along the diagonal of the simulation cubes ex- 
tends up to the redshift of interest i.e., the redshift covered 
by the first-epoch VVDS data ( z=1.5). In this simulation 
the mass-particle resolution is 2.2 • 10 12 /i~ 1 M Q and the 
present epoch is defined by a linear rms density fluctua- 
tion in a sphere of radius 8 h~ 1 Mpc of as = 0.9. 

The mass density contrasts in the redshift perturbed 
comoving coordinates S(yi,R) have been calculated at 
random positions in the simulation volume, by smooth- 
ing the particle distribution with a spherical top hat win- 
dow of length R = 8 /i^Mpc. Mass variances in different 
redshift bins are then derived using eq. i|15|) . The result 
is compared to the prediction of eq. (fTTTTf) in Fig. 00 Note 
that, even if it is clear that measurements suffer from cos- 
mic variance due to the relatively small volume sampled at 
each redshift, the predictions of eq. (|30|l are in agreement 
with the observed scaling of the linear mass variance. The 
magnitude of the correction with respect to the unper- 
turbed case is also evident; mass fluctuations recovered in 
redshift space on a 8 h~ 1 Mpc scale, in the redshift comov- 
ing coordinates, at 2=0.5(1.5) are ~ 25(35)% larger than 
in real comoving space (the correction factor is ~ 17% in 
the local Universe.) Fig. shows that this apparent en- 
hancement in the rms fluctuations results in a broadening 
of the mass PDF recovered in the redshift comoving space. 
Thus, the effect of peculiar velocities is to shrink overdense 
regions and to inflate underdense regions, enhancing the 
probability of having large density fluctuations (both pos- 
itive and negative). 

We finally compare, in various redshift intervals, the 
accuracy with which the lognormal mass PDF derived 
in the redshift comoving space (by using eq. I|30|) in I17fl 
approximates the PDF directly inferred from the Hubble 
volume simulation (see Fig. On a scale of 8 /i _1 Mpc, 
the agreement between the analytical and simulated mass 
PDFs is satisfactory at all redshifts. This holds true also 
when the mass PDFs recovered on R — 5 and 10 /i _1 Mpc 
scales are compared. 

Thus, with a good degree of confidence, we can use 
eq. H30JI to predict the PDF of mass fluctuations in red- 
shift distorted comoving coordinates (the same coordi- 
nates where the galaxy PDF is observed) and in a generic 
cosmological background. This allows us to speed up com- 
putation time and to frame the results about the biasing 
function in a generic cosmological model. 

6. Measuring Galaxy Biasing 

In this section we describe the method applied to deter- 
mine the relationship between galaxy and mass overdensi- 
ties. The galaxy ovcrdcnsity field 5 g depends in principle 
on various astrophysics! and cosmological parameters such 
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as spatial position (r), underlying matter density fluctu- 
ations (S), scale R with which the density field is recon- 
structed, cosmological time (z), galaxy colors, local gas 
temperature, non-local environment, etc. 

For the purposes of this study, we will rely on the fol- 
lowing simplifying theoretical assumptions: 

i) the efficiency of galaxy formation on a given cosmo- 
logical scale is sensitive only to the underlying mass dis- 
tribution. This means that the galaxy fluctuation field 
is in a reasonably tight one-to-one relationship with the 
underlying mass fluctuation field, and that the biasing 
scheme may be formally represented via the relationship 
8 g = b(z, 6, R)6. While such an approach represents a non- 
trivial step forward in understanding the properties of the 
biasing function b (if compared, for example, to constant 
parameterizations of the biasing relation), it is however 
evident that the biasing function could show, in principle, 
a more complex functional dependence. 

ii) The current theoretical understanding of how cluster- 
ing of DM proceeds via gravitational instability in the ex- 
panding Universe is well developed, i.e., the PDF of mass 
fluctuations of the real Universe can be safely derived via 
analytical models or N-body simulations (see discussion 
in §6) In particular, in what follows, we will consider a 
ACDM background mass distribution locally normalized 
to a s (z = 0) = 0.9. 

iii) The redshift distortions affect the densities of galaxies 
and mass in a similar way, i.e., there is no velocity bias 
between these two components, and galaxies follow the 
matter flow. 

6.1. The Method 

As described in §1, we derive the relationship between 
galaxy and mass overdensities in redshift space S g = S g (S) 
as the one-to-one transformation which maps the the- 
oretical mass PDF f(6) into the observed galaxy PDF 
g(S g ). A similar method to derive the biasing function 
has been proposed and te s ted u sing CD M simulations by 
ISigad Branchini fc Dekell ( 200(1) (see also lSzanudi fc Paa 
( 2004)). This same t echnique has been recently applied in 
different contexts bv lMarinoni fc Hudson! l|2002|) to derive 
the mass-to-light (M = M(L)) and the X-ray-to-optical 
(L x — L X (L)) fu nctions for a wide ma ss range of virializcd 
systems, and bv lOstriker et alJ l)2003h to explore the void 
phenomena in the context of hydrodynamic simulations. 

Using ea. HI2l andl3l we obtain the biasing function b(S) 
as the solution of the following differential equation 



^(-1)=-1 

b'(5)5 + b(S) = f(5)g(S g ) 



(31) 



where the prime denotes the derivative with respect to 
5, f(S) and g(S) arc the PDF of mass and galaxy fluc- 
tuations respectively, and the initial condition has been 
physically specified by requiring that galaxies cannot form 
where there is no mass. 



With this approach, we loose information on a possi- 
ble stochasticity characterizing the biasing function. The 
advantage is that we can provide a measure, on some char- 
acteristic scales R, of the local, non-linear, deterministic 
biasing function (eq. J5J) over the continuous redshift in- 
terval 0.4< 2<1.5. 

We have obtained the biasing function b(S) by numer- 
ically integrating the differential equation (|31|l . i) in dif- 
ferent redshift intervals in order to follow the evolution 
of b(S) as a function of cosmic time, and ii) using mat- 
ter and galaxy PDFs obtained by smoothing the density 
fields on R=5,8, and 10 h~ 1 Mpc in order to test the scale 
dependence of the galaxy biasing function. 

The information contained in the non-linear function 
b(S) can be compressed into a single scalar which may be 
easily compared to the constant values in term of which 
the biasing relation is usually parameterized (see eq^). 
Since, by definition, (b(d)d) = 0, the most interesting lin- 
ear bias estimators are associated to the second order mo- 
ments of the PDFs, i.e., the variance (<5^) and the covari- 
ance ( b~„ 5) . Following the prescriptions of iDekel fc Lahavl 
(1999), we characterize the biasing function as follows: 



b = 
and 



b\ 



(32) 



(33) 



where the parameter b, measuring the slope of the linear 
regression of 5 g on 5, is the natural generalization of the 
linear bias parameter defined in equation ^ and b\ is an 
"unbiased estimator" of the linear biasing parameter de- 
fined as £ 9 = fr 2 £, when the bias relation is deterministic, 
i.e., non-stochastic. The ratio r = b/bl is the relevant mea- 
sure of nonlinearity in the biasing relation; it is unity for 
linear biasing, and it is either larger or smaller than unity 
for nonlinear biasing. 

The errors in the measured values of the biasing pa- 
rameter bh have been computed using independent mock 
catalogs which implement all the selection functions of the 
VVDS. This allows us to incorporate in our error estimates 
the uncertainties due to cosmic variance. 

6.2. Testing the Method 

Before applying the biasing computation scheme (eq. I31f) 
to VVDS data, we have tested that the method can be 
meaningfully applied, i.e., it is free of systematics when 
implemented with samples of simulated galaxies which 
mimic all the observational systematics of our sample. 

The procedure consists in computing the biasing func- 
tion 5 S = b(Sp)S p between the density field S s recon- 
structed using an s-sample (representing the pseudo- 
survey sample, see §4.1) and the density fluctuations S p 
of the corresponding p-sample (representing the pseudo- 
real Universe). We have already determined the range of 
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Fig. 10. The simulated biasing function (solid-line) at differ ent cosmic epochs , between the density field traced by the 
s-sample (GALICS data simulating the VVDS sample, see IPollo et all (f2005|) and §4.1) and the density field traced 
by the p-sample (GALICS data simulating the real underlying distribution of galaxies). 6l represent the linear bias 
parameter evaluated from the nonlinear biasing function using the estimator given in eq. I|32|) . The dashed line is 
drawn at — 1 and represents the no bias case. The central cross is for reference and represents the 5 g = S = case. 
The r parameter measures the deviations from the linearity. The galaxy overdensities are reconstructed using a TH 
window of sizes R=5 /i _1 Mpc (upper panel) and R=8 /i _1 Mpc (lower panel). The shadowed area represents la errors 
in the derived biasing function. 



redshift, density contrasts and smoothing scales where 
the sample simulating all the VVDS selection functions 
(s-sample) trace the underlying density of galaxies (p- 
sample). We thus expect, for consistency, that, in that 
range, the biasing between the two samples derived by ap- 
plying our computation scheme (eq. El using the PDFs 
of the s- and p-samples) is independent of S p and equal to 
b(5 p ) = l. 

Results are presented in Fig. f° r two different TH 
smoothing scales. Note that a log-log density plot is used 
in order to emphasize the behavior of the biasing function 
in underdense regions. We conclude that on scales R>8 
h~ 1 Mpc the density recovered by a "four-passes" VVDS- 
like survey is not biased with respect to the underlying 
distribution on any density scale and in any redshift in- 
terval up to z=1.5. As a matter of fact, the linear bias 
parameter with which information contained in the bias- 
ing function can be at first order approximated is 6l ~ 
1 and the biasing relation does not show any significant 
deviation from linearity as indicated by the fact that the 
r parameter is also very close to unity. 



If the density field is smoothed on 5 /i -1 Mpc , the 
effects of the incompleteness in low-density regions (al- 
ready discussed in §4.1, see Figs. 0] and 0) become evi- 
dent. Underdense regions (log(l + 5 p ) <J —0.5) in volumes 
at redshift greater than 1 are poorly sampled with the 
VVDS survey strategy. 

In the same spirit, we have also solved eq. (1311) for 
determining the biasing relation between the PDF of the 
Hubble volume mass fluctuations and the lognormal ap- 
proximation given in eq. I|16|l . The biasing relation be- 
tween these two different descriptions of the mass density 
field is linear and consistent with the no-bias hypothe- 
sis between the two representations of the density field 
on the scales we are interested in (R > 5 /i _1 Mpc and 
log(l + S)>-l). 

7. The Biasing Function up to z ~ 1.5 

7.1. Results 

The numerical solutions of eq. I|31() for the Ai c B = — 20 + 
5 log h volume-limited VVDS sample arc plotted in vari- 
ous redshift slices, in Fig. ^2 for the cases R=8 and 10 
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1.10± 0.18 


0.97 


1.28 ±0.18 


1.4±0.4 




0.7< 2<0.9 




229 


1.12± 0.17 


0.99 


1.18 ±0.17 


1.1±0.3 




0.9< z<l.l 




289 


1.18± 0.15 


0.97 


1.17 ±0.15 


1.4±0.4 


8 


0.7< 2<0.9 


No 


1263 


0.92± 0.20 


0.97 


0.67 ±0.20 


1.6±0.4 




0.9< 2<1.1 




864 


1.03± 0.16 


0.96 


0.74 ±0.16 


1.6±0.3 




l.K 2<1.3 




440 


1.21± 0.12 


0.96 


0.82 ±0.12 


1.3±0.3 




1.3< 2<1.5 




234 


1.51± 0.10 


0.95 


0.96 ±0.11 


1.0±0.3 


8 


0.7< 2<0.9 


-18.7 


879 


0.98± 0.21 


0.98 


0.75 ±0.21 


1.6±0.4 




0.9< 2<1.1 




813 


1.01± 0.16 


0.97 


0.72 ±0.16 


1.8±0.3 


8 


0.7< 2<0.9 


-20 


279 


1.17± 0.23 


0.98 


0.89 ±0.23 


1.7±0.4 




0.9< 2<1.1 




327 


1.30± 0.17 


0.99 


0.94 ±0.17 


1.8±0.3 




l.K 2<1.3 




251 


1.33± 0.12 


0.96 


0.88 ±0.12 


1.2±0.3 




1.3< 2<1.5 




169 


1.56± 0.11 


0.95 


0.99 ±0.11 


1.1±0.3 


10 


0.7< 2<0.9 


No 


1425 


1.03± 0.22 


0.91 


0.66 ±0.22 


1.6±0.4 




0.9< 2<1.1 




955 


1.05± 0.18 


0.97 


0.64 ±0.18 


1.7±0.4 




l.K 2<1.3 




480 


1.17± 0.13 


0.90 


0.68 ±0.13 


1.3±0.3 




1.3< 2<1.5 




250 


1.55± 0.11 


0.93 


0.84 ±0.13 


1.2±0.3 


10 


0.7< 2<0.9 


-18.7 


991 


1.03± 0.25 


0.95 


0.69 ±0.25 


1.5±0.4 




0.9< 2<1.1 




900 


1.03± 0.18 


0.96 


0.63 ±0.18 


1.7±0.3 


10 


0.7< 2<0.9 


-20 


316 


1.14± 0.25 


0.92 


0.75 ±0.25 


1.6±0.4 




0.9< 2<1.1 




360 


1.26± 0.20 


0.97 


0.78 ±0.20 


1.8±0.4 




l.K 2<1.3 




266 


1.36± 0.14 


0.91 


0.78 ±0.14 


1.3±0.3 




1.3<2<1.5 




175 


1.54± 0.13 


0.93 


0.84 ±0.13 


1.3±0.3 



Table 1. Bias measurements from the VVDS first epoch data. 



/i _1 Mpc. Note that a log-log density plot is used in order 
to emphasize the behavior of the biasing function in un- 
derdense regions (note that in these units linear biasing 
appears as a curved line). 

The corresponding parameters and r (also com- 
puted for the whole flux-limited sample) are quoted in 
table ^ together with our estimates of the second mo- 
ment of the galaxy PDF an and of the skewness parameter 
S3. Both these statistics have been computed as described 
in §5. Also note that the values of a$ measured for the 
flux-limited sample are consistent with the values inde- 
pendently derived in §3.4 on the basis of the results of 
the analysis of the clustering properties of VVDS galaxies 
(Paper III). 

An empirical fit of the biasing function is obtained 
by using a fo r mula similar to the one proposed by 
iDekel fc Lahal l|l999h 



(l±a )(l±5) ai -1 <5<0 
a 2 S + a 3 S 2 S > 



(34) 



which best describes the behavior of biasing in underdense 
regions (S < 0), either the second order Taylor expansion 



of the density contrast of dark matter (Fry& Gaztanaga 
1993) 



£ 

k=0 



(35) 



which allows an easier comparison of our results with other 
studies. The best fitting parameters of these non-linear 
approximations are quoted in tabled 

The dependence of the shape of the biasing function 
on galaxy luminosity is plotted in Fig. Results are 
shown at the median depth of the VVDS sample (in 
the redshift bin 0.7< z<0.9) where faint objects (A4b < 
— 17.7 ± 51og/i) are still sampled. 

In Fig. E| we show the redshift evolution of the linear 
biasing parameter computed over the redshift interval 
0.4 < z < 1.5 for both the flux and volume limited sam- 
ples. We can conclude that biasing is not changing with 
cosmic time for z < 0.8, while there is a more pronounced 
evolution of biasing in the redshift interval [0.8,1.5]. In 
particular, the difference between the value of &l at red- 
shift z~1.5 and z~ for a population of galaxies with 
luminosity M B < -20 ± 51og/i is Ab L =~ 0.5 ± 0.14, 
thus significant at a confidence level greater than 3a. 
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R Redshift M.% a Oi a 2 a 3 b &i b 2 

h~ 1 Mpc range 



8 


0.7< z<0.9 

0. 9<z<l.l 

1. KKl.3 
1.3<z<1.5 


No 


0.36 ±0.03 
0.35 ±0.03 
0.39 ±0.04 
0.54 ±0.05 


1.36 ±0.11 
1.36 ±0.40 
1.64 ±0.13 
2.34 ±0.30 


0.60 ±0.12 
1.08 ±0.10 
1.14 ±0.11 
1.46 ±0.14 


-0.01 ±0.02 
-0.08 ±0.02 
-0.08 ±0.02 
-0.11 ± 0.02 


0.10 ±0.06 
0.23 ±0.06 
0.24 ±0.07 
0.37 ±0.12 


0.88 ±0.10 
1.18 ±0.09 
1.26 ±0.10 
1.57 ±0.15 


-0.12 ±0.08 
-0.20 ±0.06 
-0.22 ±0.06 
-0.24 ±0.08 


8 


0.7< z<0.9 

0. 9<z<l.l 

1. KKl.3 
1.3<z<1.5 


-20 


0.25 ±0.05 
0.40 ±0.06 
0.45 ±0.04 
0.48 ±0.05 


1.67 ±0.20 
2.12 ±0.33 
1.86 ±0.13 
2.51 ±0.30 


1.25 ±0.16 
1.20 ±0.16 
1.29 ±0.13 
1.46 ±0.16 


-0.09 ±0.03 
-0.05 ±0.02 
-0.11 ± 0.02 
-0.10 ±0.02 


0.18 ±0.10 
0.23 ±0.17 
0.31 ±0.08 
0.33 ±0.14 


1.30 ±0.15 
1.29 ±0.16 
1.40 ±0.13 
1.55 ±0.15 


-0.20 ±0.08 
-0.12 ±0.08 
-0.26 ±0.06 
-0.22 ±0.08 


10 


0.7< z<0.9 

0. 9< z<l.l 

1. KKl.3 
1.3<z<1.5 


No 


0.25 ±0.03 
0.18 ± 0.02 
0.39 ±0.04 
0.54 ±0.04 


0.85 ±0.10 
1.20 ± 0.08 
1.30 ±0.11 
2.09 ±0.13 


1.44 ±0.21 
1.23 ± 0.13 
1.04 ±0.13 
1.41 ±0.15 


-0.23 ±0.06 
—0.18 ± 0.04 
-0.08 ± 0.03 
-0.14 ±0.03 


0.17 ±0.07 
0.20 ± 0.04 
0.26 ±0.06 
0.39 ± 0.08 


1.20 ±0.15 
1.19 ± 0.09 
1.18 ±0.10 
1.55 ±0.12 


-0.30 ±0.12 
-0.34 ± 0.08 
-0.22 ±0.08 
-0.32 ±0.08 


10 


0.7< z<0.9 

0. 9<z<l.l 

1. KKl.3 
1.3<z<1.5 


-20 


0.21 ±0.03 
0.18 ±0.04 
0.38 ±0.04 
0.22 ±0.05 


1.33 ±0.10 
1.51 ±0.15 
1.37 ±0.17 
1.76 ±0.20 


1.45 ±0.16 
1.38 ±0.16 
1.50 ±0.15 
1.83 ±0.16 


-0.22 ±0.04 
-0.13 ±0.03 
-0.13 ±0.03 
-0.18 ±0.03 


0.18 ±0.06 
0.20 ±0.08 
0.37 ±0.10 
0.33 ±0.12 


1.26 ±0.11 
1.36 ±0.14 
1.50 ±0.14 
1.73 ±0.17 


-0.28 ±0.10 
-0.26 ±0.08 
-0.26 ±0.08 
-0.34 ±0.08 



Table 2. Best fitting parameters of the non linear biasing models given in eqs. (|r?4|t and 1)35(1 . Errors do not include 
cosmic variance. 



In Fig. ^] we show the dependence of the linear bias- 
ing parameter on galaxy luminosity. Intrinsically brighter 
galaxies are more strongly biased than less luminous ones 
at every redshift and the dependence of biasing on lumi- 
nosity at z~0.8 is in good agreement with wha t is ob- 
served in the local Universe IjNorberg et al.l 120011 . 

Given the difference in the rest-frame colors of ellipti- 
cal and irregular galaxies and the fact that the observed I 
band corresponds to bluer rest-frame bands at higher red- 
shift, the relative fraction of early- and late- type galaxies 
in our I band limited survey will change as a function of 
redshift. 

Specifically, the observed difference in th e B-band lu- 
mino sity function of early- and late-types l|Zucca et all 
120051) . implies that the VVDS survey selects preferentially 
late- type galaxies at higher redshift. It is known that at 
2=0 late-type galaxies cluster less strongly than early- 
types (e.g., Giovanelli, Haynes & Chincarini 1986, Guzzo 
et al. 1997, Giuricin et al. 2001, Madgwick et al. 2002, 
Zehavi et al. 2002), and, thus, we might observe a varia- 
tion of the amplitude of density fluctuations at high red- 
shifts just because the morphological composition of our 
sample changes. 

In order to disentangle the spurious morphological con- 
tribution to the observed evolution of the global biasing 
function we have splitted our sample according to rest 
frame colors, selecting a red ((B — I)q > 1.5; 849 galaxies 
in the 4-passes region with z > 0.7) and a blue subsample 
of galaxies ((B — I)q < 1; 1891 galaxies with z > 0.7). 
These color cuts roughly correspond to selecting, respec- 
tively, morphological types < 77 and IV according to the 
classification scheme devised by Zucca et al. 2005 for the 
VVDS sample. 

Clearly, this subsample selection does not correspond 
to the ideal case of a redshift survey sampling galaxies 
according to their rest-frame colors; however, useful in- 



formation about differences in clustering between red and 
blue populations can still be inferred. 

Note that the hypothesis on which the technique of 
comparing mass and galaxy density distributions is based 
(§6) can be straightforwardly generalized to compute 
the biasing between the density distributions of different 
galaxy types. In particular, we assume that the large scale 
velocities of late and early types are not dissimilar rela- 
tive to each other (as it is effectively observed at z=0 
e.g., iDekell Il994l iMarinoni et all Il998ll i.e., the two ve- 
locity fields are noisy versions of the same underlying field. 

Results about the color dependence of biasing are sum- 
marized in table and graphically presented in Fig. 1151 
The red sample is systematically a more biased tracer of 
mass than the blue one in every redshift interval investi- 
gated (i.e., b r > b b ), but the relative biasing between the 
two populations is nearly constant (b r /b b ~ 1.4 ± 0.1) 

7.2. Analysis and Discussion 

7.2.1. Biasing for the Global Galaxy Population 

Here we examine and interpret the results derived in 
the previous section. We begin by discussing the general 
shape of the non-linear biasing function, for the global 
population, in different density regions. Our results can 
be summarized as follows: 

i) in underdense regions (1 + S <1) the local slope of 
the biasing function b(S) is always larger than unity even 
when the global slope is b^ < 1 (see for example Fig. 
I12fl . The fact that galaxies in low-mass density regions 
are always positively biased with respect to the mass dis- 
tribution (i.e., locally b > 1) is possibly physically caused 
by the fact that galaxies do not form in very low-density 
mass regions, i.e., below some finite mass underdensity 
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Fig. 11. The observed biasing function (solid-line) recovered for the density field smoothed on scales R= 8 ft, _1 Mpc 
(upper panel) and 10 h~ 1 Mpc (lower panel) and for different redshift bins (from left to right) in the volume-limited 
VVDS sample (M.% — — 20 + 51og/i). The dotted line represents the linear biasing model 8 g = 6l<5 while the no-bias 
case (bL = 1) is shown with a dashed line. The central cross is for reference and represents the S g = S = case. The 
shaded area represents la errors in the derived biasing function. Errors take into account the noise in the observed 
galaxy PDF (g(S g )), but do not include uncertainties due to cosmic variance. 



M-51og h<-17.7 M-51ogh<-18.7 M-51og h<-20 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

log(1 +6 ) log(1 +6 ) log(1 +6 ) 



Fig. 12. The biasing function (solid-line) on scales R= 8 h 1 Mpc and in the redshift interval 0.7< z<0.9 computed 
for different luminosity classes. Symbols are as in Fig. 1111 



the galaxy formation efficiency drops to zero. Using the 
biasing relation given in eq. 134(1 the characteristic mass 
density threshold S c below which very few galaxies form 
(5 g < —0.9), can be approximated as 

i n , x \ l + log(l + ao) 

log(l + S c ) (36) 



There is evidence that this mass-density threshold, 
characterizing regions avoided by galaxies, increases as 
a function of redshift (see Fig. Illfl and luminosity (see 
Fig. 1121) . If we consider R=8 /i _1 Mpc and the redshift bin 
0.7< z <0.9 we see that while faint galaxies seem to be 
present even where the mass density contrast is very low 
(left panel of Fig. El log(l + S c ) = -0.96 ± 0.09), brighter 
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Linear bias evolution 
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Fig. 13. The redshift evolution of the linear biasing pa- 
rameter 6l for the volume- limited (A4 C B <-20+51ogh) sub- 
sample (filled squares) is compared to the evolution of the 
biasing parameter for the whole flux-limited VVDS-02h 
sample (empty squares). Since there is no significant evi- 
dence of scale dependence in the biasing relation, we have 
averaged the biasing parameters measured on 5,8, and 10 
h~ 1 Mpc scales in order to cover the full redshift baseline 
0.4< 2<1.5. For clarity, only the errorbars corresponding 
to the volume-limited sample are shown. The triangle rep- 
resents the z~0 bias inferred for 2dFGRS galaxies hav- 
ing median L/L* ~ 2 (i.e., the median luminosity of the 
volume-limited VVDS sample) as explained in th text. 



N 

n 




Fig. 14. Comparison between the galaxy linear bias pa- 
rameter measured in the redshift interval 0.4 < z < 0.9 
for 3 different luminosity classes (squares) and the cor- 
responding local estimates provided by the 2dFGRS (tri- 
angles). Points with increasing sizes correspond to three 
different volume-limited VVDS subsamples, i.e Mb — 
51ogft < -17.7, < -18.7 and < -20, respectively. For 
clarity, squares with increasing size have been progres- 
sively displaced rightward to avoid crowding. The z~0 
measurements have been interpolated by using the for- 
mula describing the luminosity depen dence of the 2dFGRS 
bias parameter l|Norberg et alll200l[) . the bi as parameter 
for th e 2dFGRS L* sample (i.e., 6* = 0.92 l|Verde et all 
2002)) and the median luminosity of the three VVDS sub- 
samples (L/L*=0.52,0.82,2.0 respectively). 



galaxies do not seem to form in deep mass underdensi- 
ties (right panel of Fig. El log(l + 8 C ) = -0.73 ± 0.11). 
Therefore low-density regions are preferentially inhabited 
by low luminosity galaxies. 

Moreover the mass-density threshold below which the 
formation of bright galaxies (A4b < — 20 + 5 log ft.) seems 
to be inhibited increases, irrespective of the scale investi- 
gated (see Fig. IIH as a function of redshift. On a scale 
R=8 /i _1 Mpc, the threshold shifts from log(l + 5 C ) = 

-0.73 ± 0.10 at z=0.8 to log(l + S c ) 0.55 ± 0.07 at 

z=1.4 This suggests that galaxies of a given luminosity 
were tracing systematically higher mass overdensities in 
the early Universe, i.e, as time progresses, galaxy forma- 
tion begins to take place also in lower density peaks. 

ii) Even in regions where the mass density distribution 
is close to its mean value (l+<5 ~ 1) bright galaxies are 
not unbiased tracer of the mass-overdensity field (Fig. II 1(1 . 
This can also be seen by setting 5 = in ea.l3*4land noting 



that S g (S = 0) = ao > for both analyzed samples (flux- 
and volume-limited) in all redshift ranges (see table 2). 
This result is at variance with what is expected within 
the simple linear biasing picture, where, by construction, 
5 g (8 = 0) = 0. 

iii) In higher matter-density environments (l+<5 > 1) 
galaxies were progressively more biased mass tracers in 
the past, i.e., the local slope b(S) systematically increases 
with redshift on every scale investigated (Fig. lllfl . There is 
some indication that, at the upper tail of the mass density 
distribution, galaxies are anti-biased with respect to mass 
on all scales (i.e., the local slope is b(5) < 1 for 8 ^S> 1). 
Antibiasing in overdense regi mes is a feature actually 
observed in simulations (e.g. , ISigad. Branchini fc D ekel , 
l2000HSomerville et all l200llj and e xpected in theoretical 
models (e.g., iTaruva fc Satll2000j) . Physically this could 
be due to the merging of galaxies which reduces the num- 
ber density of visible objects in high density regions or 
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because galaxy formation is inhibited in regions where the 
gas is too hot to collapse and form stars. 

iv) In general the linear approximation offers a poor 
description of the richness of details encoded in the biasing 
function. As a matter of fact the linear biasing function 
(dotted line in Fig. Illland ll2|) poorly describes, in many 
cases, the observed scaling of the biasing relation (solid 
line). At the comoving scales of R=5, 8 and 10 /i~ 1 Mpc, 
non-linearities in the biasing relation are typically (1— r) <J 
10% in the redshift range investigated. We find that the 
ratio 62/61 between the quadratic and linear term of the 
series approximation given in eq. I)35|l is nearly constant 
in the redshift range 0.7 < z < 1.5 and does not depend 
on luminosity (i.e. it is nearly the same for the flux- and 
volume-limited subsamples) or smoothing scale. We find 

that, on average, b 2 /b 1 0.15 ± 0.04 for R = 8 /i _1 Mpc 

and b 2 /b 1 0.19 ± 0.04 for R = 10 h^Mpc. . 

To facilitate comparison with other studies, which gen- 
erally focus on the linear representation of biasing, we now 
discuss the properties of the linear approximation of our 
biasing function. The general characteristics of the linear 
parameter bh can be summarized as follows: 

v) by inspecting tabled we do not find any significant 
evidence that the global value of the linear biasing param- 
eter bh depends on the smoothing scale. Any possible sys- 
tematic variation, if present, is smaller than the amplitude 
of our errorbars (~0.15). This scale independence in the 
biasing relation extends into the high redshift regimes sim- 
ilar conclusions obtained in the local Universe b y the 2dF- 
GRS on scales >5 /i^Mpc l|Verde et all 120021) . Moreover 
our results may be interpreted as a supporting evidence for 
theoretical arguments suggesting that bias is expected to 
be scale-independent on scales larger than a few ft,~ 1 Mpc 
(e.g., Mann, Peacock & Heavens 1998, Weinberg et al. 
2004). 

Since we find no evidence of scale-dependent bias, and 
since with different R scales we are probing different red- 
shift regimes, in Fig. 1131 we have averaged the linear bi- 
asing parameters measured on 5,8, and 10 ft, _1 Mpc scales 
(values quoted in Table 1) in order to follow, in a continu- 
ous way, the redshift evolution of the linear galaxy biasing 
over the larger redshift baseline 0.4< z<1.5. Fig.ll3lshows 
that &l for galaxies brighter than M.b = —20 + 51og/i 
changes from 1.10 ± 0.18 at z ~ 0.55 to 1.55 ± 0.12 at 

1.4. 

An even steeper variation is observed for the biasing 
of the flux-limited sample, indicating that biasing depends 
on galaxy luminosity. Fig. 1 131 shows that the ratio between 
the amplitude of galaxy fluctuations and the underlying 
mass fluctuations declines with cosmic time. This scal- 
ing is effectively predicted withi n the f r amew ork of the 
peaks-biasing theoretical model <|Kaiserl 119841) . At early 
times, galaxies are expected to form at the highest peaks 
of the density field since one needs a dense enough clump 
of baryons in order to start forming stars. Such high-er 
peaks are highly biased tracers of the underlying mass 
density field. According to this picture, as time progresses 
and the density field evolves, galaxy formation moves to 



lowcr-er peaks, nonlinear peaks become less rare events and 
thus galaxies become less biased tracers of the mass den- 
sity field. Additional "debiasing" mechanisms may con- 
tribute to the observed scaling shown in Fig. El It is likely 
that the densest regions stop forming new galaxies because 
their gas becomes too hot, canno t cool efficiently, and thus 
cannot collapse and form stars llBlanton et al.lll999h . As 
galaxy formation moves out of the hottest (and rarest) 
regions of the Universe, the biasing decreases. Finally, we 
also note that in order to derive the biasing function we 
have assumed that there is no difference in the veloc- 
ity field of the luminous and matter components. After 
galaxies form, they are subject to the same gravitational 
forces as the dark matter, and thus they tend to trace the 
da rk matter distr i bution more closel y with time as shown 
bv iDekel fc Reesl <|l987h :lFrl l|l99filk iTegmark fc Peebles! 
<l!998l) . 

vi) In Fig. ll3l we also show, for comparison, the value of 
the 2dFGRS linear biasing parameter inferred at z=0.17 
(the effective depth of the survey) as the ratio between 
the as value measured by Croton et al. 2004 (in redshift- 
distorted space; see their Fig 3 and 4) for a sample of 
objects with —21 < A4b — 51og/i < —20 (which actu- 
ally brackets the median luminosity of our volume lim- 
ited sample ~ 2L*), and the rms of mass fluctuations 
(in redshift-distorted space) in a ACDM background (see 
§6). This value (~ 1.07 ± 0.06) is in excellent agreement 
with what one would independently o btain by combining 
the linear bias parameter measured bv I Verde" et all <|2002l) 
for the whole 2dFGRS (1.04 ±0.11 ) with the bias scaling 
law recipe of iNorberg et al ] l|200ll) . i.e., b(z = 0.17, L = 
2L*) = 1.07 ±0.13. 

We can conclude that the time dependence of biasing is 
marginal (db/b ~ 7±25%) for z<0.8 while it is substantial 
(db/b ~ 33±18%) in the resdhift interval [0.8-1.5]. The ob- 
served time evolution of bias is well described by the sim- 
ple scaling relationship b L = 1 + (0.03±0.01)(1 + z) 3 ' 3±0 6 
in the interval 0< z <1.5. 

Assuming a linear biasing scheme, one may note that 
this result was already implicit in Fig. of §4. The rms 
fluctuations of the mass density field on a 8 h~ 1 Mpc scale 
decrease monotonically with redshift by a factor of ~ 22% 
and ~ 23% in the redshift intervals [0.17-0.8] and [0.8-1.4], 
respectively; thus, a nearly constant bias is predicted in 
the redshift range z=[0.17-0.8] because the rms fluctua- 
tions of the galaxy density field are also decreasing by a 
factor ~ 16% in this same interval. Since, instead, as of 
galaxies is marginally increasing in the range z=0.8-1.4 
(da I a ~ 10%, see table 1), over this redshift baseline the 
biasing evolves rapidly. 

vii) Bright galaxies are more biased mass tracers than 
the general population (see Fig. I12|l . This result confirms 
and extends into the high redshift domain the luminosity 
dependence of biasing which is observed in local samples 
of galaxies (e.g., Benoist et al. 1996, Giuricin et al 2001, 
Norberg et al. 2001, Zehavi et al. 2002). Specifically, in 
Fig. 1141 we show the dependence of galaxy biasing from 
luminosity measured in the redshift interval 0.4<z<0.9 
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using three different volume-limited VVDS subsamples 
(i.e., M B - 5\ogh < -17.7, < -18.7 and < -20 re- 
spectively) and compare their linear biasing parameters 
with those observed locally for a sample of objects hav- 
ing the same median luminosities of the VVDS subsam- 
ples (i.e., L/L*=Q. 52,0.82,2.0 respectively). The local esti- 
mates have been computed on the basis of t he scaling rela- 
tionsh ip b/b*= 0.85 + 0.15 L/L* derived bv lNorberg et all 
using the 2dFGRS sa mple, assuming the b* value 
given by IVerde et alJ l|2002|) . As shown above for the 
volume-limited sample, no significant evolution is seen up 
to z^0.8 also when the dependence of bias from luminos- 
ity is analyzed. 

Finally, we note that, as already discussed in §5, galax- 
ies with the same luminosity at different redshifts may 
actually correspond to different populations. Since, as we 
have shown, biasing increases with luminosity also at high 
redshift, and since the measured value of M* for our sam- 
ple at redshift z=0.4(1.5) (Paper II) is fainter (brighter) 
than the cut-off magnitude M C B — — 20 + 51og/i, we can 
infer that for a population of objects selected, at 

any given redshift, in a narrow luminosity range around 
A4*(z) should increase with redshift even more than what 
we have measured for our volume-limited sample (sec Fig. 
113(1 . A more detailed analysis of the biasing for M*(z) 
galaxies will be presented in the future, when a larger 
VVDS data sample will be available. 



7.2.2. Biasing as a Function of Galaxy Color 

Results summarized in table |3 and presented in Fig. [TBI 
show that, on scales R=8 ft. _1 Mpc, the red sample is a 
more biased tracer of mass than the blue one in every 
redshift interval. Similarly to what we have found for the 
global population, there is some indication of a systematic 
increase as a function of redshift of the biasing of bright 
red and blue objects even if, because of the large errorbars, 
this trend is not statistically significant. 

We can compare our results to the biasing measured 
for extremely red objects (EROS), i.e., objects with ex- 
tremely red colors ((R — K) vega > 5). Using the results of 
the correlation analysis of Firth et al. (2002), we obtain, 
for their (/ — H) vega > 3, H vega < 20.5 sample (which 
has a median blue luminosity Mb = —20.3 + 51og/i), 
bf ROS {z ~ 1.2) - 2.3 ± 0.6. Considering the results of 
Daddi et al. (2001), who analyzed a sample of EROS 
with (R — K) vega > 5 (which roughly corresponds to 
(/ — H) vega > 3), K vega < 19.2 sample, we conclude that 
bf ROS ( z ~ 1.2) ~ 4 ± 1. These values for the galaxy bi- 
asing are respectively ~ 0.5 and 1.8cr higher than that 
measured for our sample of bright (Mb < — 20 + 51og/i) 
but moderately red galaxies (b r L (z ~ 1.2) = 2 ± 0.5). One 
may interpret this results as an indication for the reddest 
objects being more strongly biased then moderately red 
galaxies of similar luminsity. Anyway, given the large er- 
rorbars, the evidence that, at z ~ 1.2, the biasing proper- 
ties of these two differently selected populations are differ- 



Redshift 


Volume 


b rel 


bl 


bf 


range 


limited 








0.7<z<0.9 


No 


1.3 ±0.2 


1.3 ±0.5 


1.0 ±0.5 


0.9<z<l.l 


No 


1.3 ± 0.2 


1.4 ±0.5 


1.1 ± 0.4 


l.Kz<1.3 


No 


1.4 ±0.2 


1.6 ±0.5 


1.2 ± 0.4 


1.3<z<1.5 


No 


1.4 ±0.3 


2.3 ± 0.4 


1.6 ±0.4 


0.7<z<0.9 


-20 


1.4 ±0.3 


1.5 ±0.6 


1.1 ±0.6 


0.9<z<l.l 


-20 


1.3 ±0.3 


1.6 ±0.6 


1.2 ± 0.5 


l.Kz<1.3 


-20 


1.5 ± 0.2 


2.0 ±0.5 


1.3 ± 0.4 


1.3<z<1.5 


-20 


1.3 ±0.3 


2.1 ± 0.4 


1.6 ± 0.4 



Table 3. The biasing parameters for red and blue VVDS 
subsamples on a scale R=8 /i _1 Mpc. 



ent is not statistically significant. As a matter of fact, the 
values quoted above are also consistent with an alterna- 
tive hypothesis, i.e., the strength of the EROS fluctuations 
with respect to the mass fluctuations is not exceptional 
when compared to the density fluctuations observed in a 
sample of high redshift, moderately red galaxies of similar 
luminosity. 

The specific values of the biasing parameter at each 
cosmic epoch are affected by large errors due to the sparse- 
ness of our volume-limited subsamples, and to the pres- 
ence of cosmic variance. One way to bypass uncertainties 
due to cosmic variance consists in computing the rela- 
tive biasing function b rel (8) — b r (5)/b b (5) between the red 
and blue subsamples. As the subsamples are drawn from 
the same volume, this ratio should be minimally affected 
by the finiteness of the volume probed by the first epoch 
VVDS data. 

Results about the relative biasing between galaxy of 
different colors are graphically shown in the lower panel 
of Fig. El while estimates of the corresponding bh are 
quoted in table [21 

We do not observe any trend in the relative biasing 
between red and blue volume-limited subsamples in the 
redshift range 0.7<z<1.5. Moreover, our best estimate 
Wf 1 ~ 1.4±0.1 is in excellent agreement with what is found 
for nearly the same color-selected populations both locally 
(Willmer et al. (1998) found that, on a scale R=8 /i _1 Mpc, 

b r L el = b({B - R) , vega > 1.3)/b({B - R) ,vega < 1.3) = 

1.4 ± 0.3, while Wild et al. 2005 using the 2dFGRS found 
on the same scale b r L el = b((B - R) , vega > 1. 07)/b((B- 
R)n , ? „, < 1.07) = 1.5 ± 0.07) and at z~l ijCoil et alJ 
(2004) found, on a scale R=8 h^Mpc, that b r L el = b{{B - 
R)o > 0.7)/b((B-R)o < 0.7) = 1.41±0.10). Thus, VVDS 
results suggest that there is no-redshift dependence for 
the relative biasing between red and blue objects up to 
z ~ 1.5. Possible systematics could conspire to produce 
the observed results; the linear approximation may not 
always captures, in an accurate way, all the information 
contained in the biasing function, and more importantly, a 
purely magnitude limited survey samples the red and blue 
populations at high redshift with a different efficiency (see 
discussion in §7.1). 



24 



Marinoni et al.: The VVDS: Galaxy biasing up to z=1.5 



Smoothing R = 8h"'Mpc 



2.5 h 
2.0 
1.5 
1.0 h 
0.5 '- 



0.0 
1.81- 
1.6 
1.4 
1.2 
1.01- 
0.8 t 



(B-l)>1.5 vol-limited M<-20 + 5log h 
(B — 1)< 1 .0 vol-limited M<-20 + 5log h 

—M 1 ■ ■ ■ 1 ■ . ■ 1 ■ ■ ■ 1 ' 



x 



0.6 



0.8 



1.0 



1.2 



1.4 



1.6 



Fig. 15. Upper panel: Redshift evolution of the galaxy 
bias on a scale R=8 h~ 1 Mpc for the red (squares) and 
blue (triangles) galaxies in the volume limited samples. 
For clarity, the triangles have been slightly displaced right- 
ward to avoid crowding. Black diamonds represent the 
global bias for galaxies brighter than M C B = —20 4- 5 log h. 
Lower panel: the relative bias between the red and blue 
population (b rel = V /b b ) is shown as a function of red- 
shift. The filled and shaded areas represents the la confi- 
dence region of the z~0 value for the relative bias derived 
by Wilmer et al. (1998) and Wild et al. (2005) respec- 
tively. The diamond represents the relative bias measured 
by Coil et al. (2004) in the redshift interval (0.7-1.35). 



In principle, the relative bias could be further studied 
as a function of scale. For example, locally, there is ev- 
idence of scale dependence in the relative bias with the 
bias decreasing as scale increases (Willmer et al. 1998, 
Madgwick et al. 2003, Wild et al. 2005). However, the 
sample currently available is not sufficiently large to ob- 
tain proper statistics on this effect, although this should 
be measurable from the final data set. 

Finally, we note that no differences in the value of 6£ ei 
are seen by comparing volume-limited subsamples with 
the flux-limited one in different redshift intervals (see table 

EJ). 

Thus, we can deduce that in each redshift bins 
b R {M B < -20)/b R ~ b B (M B < -20)/b B . In other terms 
the biasing between the most luminous objects of a par- 
ticular color and the global population of objects of the 
same type appears to be independent of galaxy colors (see 
table 0). 



8. Comparison with Theoretical Predictions 

In this section we compare our results about the biasing 
of the A4 B = —20 + 51og/i volume-limited, global galaxy 
sample, with predictions of different theoretical models. 

Since we have found that the distribution of galaxy and 
mass fluctuations are different and the bias was system- 
atically stronger in the past, we can immediately exclude 
the scenario in which galaxies trace the mass at all cosmic 
epochs. We thus consider more complex theoretical de- 
scriptions of the biasing functions, in particular three dif- 
ferent pictures based on orthogonal ideas of how evolution 
proceeds: the conserving, th e merging, and the star f orm- 
ing biasing models (see e.g.. lMoscardini et alJ l|l998() \ 

In the first model th e number of galaxies is conserve d 
as a function of time l|Dekel fc Reesl Il987t IfTvL Il996ft . 
This model does not assume anything about the distri- 
bution and mass of dark matter halos or their connection 
with galaxies. In this scheme one assumes that galaxies 
are biased at birth and then they follow the flow of mat- 
ter without merging, in other terms they behave as test 
particles dragged around by the surrounding density fluc- 
tuations. Because the acceleration on galaxies is the same 
as that on the dark matter, the gravitational evolution af- 
ter formation will tend to bring the bias closer to unity, as 
described by Fry (1996) and Tegmark & Peebles (1998). 

T he evolution of the bias is given by 
(e.g., iTeemark fc Peeblesl fl998^ 



(37) 



where 6/ is the bias at the formation time Zf. 

An alternative picture for the bias evolution, which 
explicitly t akes into account g alaxy merging, has been 
proposed bv lMo fc White! i|199(tI) who gave analytical pre- 
scriptions for computing the bias of halos using the Press 
& Schechter formalism. 

If we explicitly assume that galaxies can be identified 
with dark matter halos, an approximate expression for the 
biasing of all halos of mass > M existing at redshift z (but 
which collapsed at redshift greater than the observation 
redshift, see discussion in Matarrese et al. 1997) is given 

by 



b(M, z) 



1 

1 + ^r- 



5 c \a 2 (M,z) 



1 



(38) 



where S c ~ 1.69 is the linear overdensity of a sphere which 
collapses in an Einstein-de Sitter Universe and c(M, z) is 
the linear rms fluctuations on scales corresponding to mass 
M at the redshift of observation. 

The third model is also framed within the peaks- 
biasing formalism. It assumes that the distribution of 
galaxies with luminosity > L is well traced by halos with 
mass >M, and predicts the biasing of objects that just 
collapsed at the redshift of observation (e.g., Blanton et 
al. 2000). In this star forming model, 



b(M, z) = l + 



a 2 (M,z) 



(39) 
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Fig. 16. The redshift evolution of the linear biasing pa- 
rameter bh for the volume-limited (M% = — 20 + 5\ogh) 
sample (see Fig 13) is compared to various theoretical 
models of biasing evolution. The dotted line indicates 
the conserving model normalized at b/(z = 1.4) = 1.28, 
the solid and dashed lines represent the star forming and 
merging models with the mass thresholds set at 3.2 • 10 11 
and 2.4 • 10 12 /i _1 M Q respectively. 



Model 


Best 


fitting parameters 


X 2 /dof 


Conserving 


b f (z 


= 1.4) = 1.28 ±0.03 


2 


Merging 


M = 


2.4 ■ 1O 12 /i _1 M 


5.5 


Star forming 


M = 


(3.2 ±3) ■ 10 10 /i _1 M© 


0.7 



Table 4. Best fitting parameters and the corresponding 
X 2 values for various biasing models 



represents the biasing of galaxies that formed in a nar- 
row time interval around redshift z (i.e., galaxies which 
experienced recent star formation at redshift z.) 

Clearly the above models, are based on a set of theo- 
retical ingredients which represent a crude approximation 
of the complex multiplicity of physical phenomena enter- 
ing the cosmic recipe of galaxy biasing. In this context, 
our goal is to investigate the robustness of the simplifying 
assumptions on which theoretical models are based, and 
explore the validity or limits of their underlying physical 
motivations. 

Theoretical predictions are compared to observations 
(VVDS data plus the local normalization derived from 
2dFGRS data) in Fig. The best fitting parameters 
for each model are evaluated using a x 2 statistics and 
are quoted, together with the corresponding minimum x 2 
value of the fit, in table 4. 

The best fitting galaxy conserving model is obtained 
when the bias at birth is bt(zf = 1.4) = 1.28±0.03 and the 



corresponding normalized x 2_va hie is x% — 2. As shown 
in Fig. the redshift evolution predicted by this model is 
much weaker than suggested by data. Thus, the gravita- 
tional debiasing is a physical mechanism that alone may 
not fully explain the observed redshift evolution of the bi- 
asing, in the sense that it significantly underpredicts the 
rate of evolution. 

The redshift evolution is more pronounced in the 
merging model (specifically, in Fig. 1161 we show the 
bias evolution of galaxies hosted in halos having mass 
M £ 2.4 • 10 12 /i _1 A/ Q ). While this model successfully 
des cribes the redshift dependence of the bias ing of ha- 
los l|Mo fc Whitd . Il996l ISomerville etafl l200l|) it poorly 
accounts for the redshift evolution of the bias of galax- 
ies (with M.b < —20 + 51og/i) between z=0 and 2=1.5 
which is slower than predicted (x% = 5.5 for the best 
fitting model). Thus, although merging is an important 
mechanism for describing the evolution of matter cluster- 
ing, our result implies that merging processes affect galax- 
ies in a less dramatic way than halos. Since in the Press & 
Schechter formalism halos are required to merge instan- 
taneously in bigger units at the redshift of observation, 
our result would imply, also, that the merger time-scales 
of galaxies is different from that of halos. Moreover, se- 
lecting galaxies with a fixed luminosity threshold may not 
correspond, over such a wide z range as that investigated 
here, to selecting halos above a given fixed mass threshold. 
In this sense our result would be suggestive of evolution 
in the mass-to-light ratio as a function of time. 

In Fig El we also show the expected redshift evolu- 
tion for the star forming model for halos of M ^ 3.2 • 
lO lo /i _1 M0. In this case, the agreement between model 
and observations is better {x 2 N = 0.7). Clearly this does 
not mean that we are analyzing a sample of objects that 
just collapsed and formed stars at the time they were 
observed; as a matter of fact the model cannot capture 
all the physical processes shaping the biasing relation. 
Moreover, the low value fitted for the mass threshold is 
somewhat unrealistic for the bright objects we are consid- 
ering. Notwithstanding, Blanton et al. 2000 already noted 
that the prediction of this biasing model is not much dif- 
ferent from the biasing evolution expected for the general 
population of galaxies in a hydrodynamical simulation of 
the large scale structure. 

Our analysis seems to suggest the apparent need of 
more complex biasing models that better approximate 
the observed biasing evolution. Understanding our re- 
sults completely, however, will require more discrimina- 
tory power in the data, and, thus, a larger VVDS sample. 

9. Summary and Conclusions 

Deep surveys of the Universe provide the basic ingredients 
needed to compute the probability distribution function 
of galaxy fluctuations and to constrain its evolution with 
cosmic time. The evolution of the galaxy PDF may shed 
light onto the general assumption that structures grows 
via gravitational collapse of density fluctuations that are 
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small at early time. When this statistics is combined with 
analytical CDM predictions for the PDF of mass, useful 
insights into the biasing function relating mass and galaxy 
distributions can be obtained. 

In this paper, we have explored the potentiality of this 
approach by analyzing the first-epoch data of the VVDS 
survey. This is the largest, purely flux-limited sample of 
spectroscopically measured galaxies, currently available in 
a continuously connected volume and with a robust sam- 
pling up to redshifts z^l.5. It is worth to emphasize that 
the VVDS is probing the high redshift domain at I<24 
in the VVDS-02h-4 field with the same sampling rate of 
pioneer surveys of the local Universe such as the CFA (at 
z^O) and, more recently, the 2dFGRS (at z~0.1). 

Particular attention has been paid to assess the com- 
pleteness of the VVDS sample and to test the statistical 
reliability of the PDF of VVDS galaxy fluctuations. In 
particular: 

a) by applying the VVDS observational selection func- 
tions to GALICS semi-analytical galaxy simulations we 
have explored the region of the parameter space where 
the PDF of VVDS-likc densities traces in a statistically 
unbiased way the parent underlying PDF of the real dis- 
tribution of galaxy overdensities. 

b) we have reconstructed the VVDS galaxy density 
field on different scales R= [5,8,10] /i _1 Mpc by correcting 
the density estimator for various VVDS selection func- 
tions. The final density map for the flux-limited sample 
has been Wiener filtered in order to minimize the shot- 
noise contribution 

By studying the PDF of galaxies in the high redshift 
Universe we have found that the peak of the galaxy PDF 
systematically shifts to lower density contrasts as a func- 
tion of redshift and that the probability of observing un- 
derdense regions is greater at 2^0.7 than it was at z^-d.5. 
Both these effects provide strong supporting evidence for 
the standard assumption that the large scale structure is 
the result of the gravitational growth of small primordial 
density fluctuations in an expanding universe. 

First, within the paradigm of gravitational instabil- 
ity, the assembling process of the large-scale structures is 
thought to be regulated by the interplay of two compet- 
ing effects: the tendency of local self-gravity to make over- 
dense regions collapse and the opposite tendency of global 
cosmological expansion to move them apart. A key signa- 
ture of gravitational evolution of density fluctuations in an 
expanding Universe is that underdense regions, experienc- 
ing the cosmological matter outflow, occupy a larger vol- 
ume fraction at present epoch than in the early Universe. 
Secondly, both this effects, the peak shift and the develop- 
ment of a low density tail, could indicate the existence of a 
time-evolving biasing between matter and galaxies, since 
galaxy biasing, systematically increasing with redshift, of- 
fers a natural mechanism to re-map the galaxy PDF into 
progressively higher intervals of density contrasts. 

This last interpretation is confirmed by our measure- 
ments of the evolution properties of the second and third 
moments of the galaxy PDF. We find that i) the rms 



amplitude of the fluctuations of bright VVDS galaxies is 
with good approximation constant over the full redshift 
baseline investigated. Specifically we have shown that, 
in redshift space, <j$ for galaxies brighter than M. C B — 
—20 + 51og/i has a mean value of 0.94±0.07 in the red- 
shift interval 0.7< z<1.5;ii) the third moment of the PDF, 
i.e., the skewness, increases with cosmic time. Its value at 
z ~ 1.5 is nearly 2a lower than measured locally by the 
2dFGRS. Both these results, when compared to predic- 
tions of linear and second order perturbation theory, un- 
ambiguously indicate that galaxy biasing is an increasing 
function of redshift. 

Exploiting the sensitivity of the galaxy PDF to the 
specific form of the mass-galaxy mapping, we have derived 
the redshift-, density-, and scale-dependent biasing func- 
tion b(z,S,R) between galaxy and matter fluctuations in 
a ACDM Universe, by analyzing the Jacobian transforma- 
tion between their respective PDFs. Particular attention 
has been paid to devise an optimal strategy so that the 
comparison of the PDFs of mass and galaxies can be car- 
ried out in an objective and accurate way. Specifically, we 
have corrected the lognormal approximation, which de- 
scribes the mass density PDF, in order to take into ac- 
count redshift distortions induced by galaxy peculiar ve- 
locities at early cosmic epochs where the mapping between 
redshifts and comoving positions is not linear. In this way, 
theoretical predictions can be directly compared to obser- 
vational quantities derived in redshift space. 

Without a-priori parameterizing the form of the bi- 
asing function, we have shown its general non trivial 
shape, and studied its evolution as a function of cosmic 
epochs. Our main results about biasing in the high red- 
shift Universe can be summarized as follows: 
i) we detect non-linear effects in the biasing relation. The 
ratio between the quadratic and linear term of the biasing 
expansion (cfr eq. Iri5|) is different from zero at a confi- 
dence level greater than 3<r in all the redshift bins and 
for all the smoothing scales probed. This result confirms 
a general prediction of C DM-based hierarchical models of 
galaxy formation (e.g. , ISiead. Branchini fc Dekell . feoOO: 
ISomerville et allhoOlJ) . Such non- linear distortions of the 
biasing function are not o bserved locally in the 2dFGRS 
sample l|Verde et all 120021 but see Baugh et al. 2004), al- 
though indirect evidences of a non-linear bias at z ~ 
exist (Benoist et al. 1999, Baugh et al. 2004). 
it) The biasing function rises sharply in underdense re- 
gions (the local slope is b(S) > 1) indicating that below 
some finite mass density threshold the formation efficiency 
of galaxies brighter than Mb < —20 + 5 log h drops to 
zero. This threshold shifts towards higher values of the 
mass density field as the luminosity or the redshift of the 
galaxy population increases. 

in) We do not observe the imprints of scale-dependency 
in the biasing function a behavior in agree ment with re- 
sults derived from more local surveys at z~ l|Verde et all 
I2002J) . 

iv) By representing the biasing function in linear approx- 
imation, we have found that the linear biasing param- 
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eter 6^ evolves with cosmic time: it appears that we 
live in a special epoch in which the galaxy distribution 
traces the underlying mass distribution on large scales 
(b L ~ 1), while, in the past, the two fields were progres- 
sively dissimilar and the relative biasing systematically 
higher. The difference between the value of b^ at redshift 
z^l.5 and z~ for a population of galaxies with luminos- 
ity Mb < —20 + 5 log ft, is significant at a confidence level 
greater than 3er (Afo^ ~ 0.5 ± 0.14). In this interval, the 
essential characteristics of the time evolution of the linear 
bias are well described in terms of the phenomenological 
relationship b L = 1 + (0.03 ± 0.01)(1 + z ) 3 - 3+0 - 6 . 

v) Over the redshift baseline investigated, the rate of bi- 
asing evolution is a function of redshift: z~0.8 is the 
characteristic redshift which marks the transition from 
a "minimum-evolution" late epoch to an early period 
where the biasing evolution for a population of M.b < 
—20 + 5 log h galaxies is substantial (~ 33 ± 18% between 
redshift 0.8 and 1.5). 

vi) Brighter galaxies are more strongly biased than less 
luminous ones at every redshift and the dependence of 
biasing on luminosity at z^0.8 is in good agreement with 
what is observed in the local Universe. 

vii) By comparing our results to predictions of theoreti- 
cal models for the biasing ev olution, w e have shown that 
the galaxy con serving model (|Frvl fl996) and halo merging 
llMo fc Whitelll99d) ' model offer a poor description of our 
data. This result could suggest that the gravitational debi- 
asing and the hierarchical merging of halos may not be the 
only physical mechanisms driving the evolution of galaxy 
biasing across cosmic epochs. At variance with these re- 
sults, the star forming model (Blanton et al. 2000) seems 
to describe better the observed redshift evolution of the 
linear biasing factor. 

viii) After splitting the first-epoch redshift catalog into 
red and blue volume-limited subsamples, we have found 
that the red sample is systematically a more biased tracer 
of mass than the blue one in every redshift interval inves- 
tigated, but the relative biasing between the two popula- 
tions is nearly constant in the redshift range 0.7< z<1.5 
(b r /b b ~ 1.4 ± 0.1), and comparable with local estimates. 
Moreover, we have found that the bright red subsample is 
biased with respect to the general red population in the 
same way as the bright sample of blue objects is biased 
with respect to the global blue population thus indicating, 
that biasing as a function of luminosity might be, at first 
order, independent of colors. 

ix) Because the VVDS and various EROS samples are not 
yet large enough, the bias of our sample of bright and 
moderately red objects at z~l is not statistically dissim- 
ilar from that expected for EROS of similar luminosity, 
even if the EROS biasing appears to be systematically 
larger. 

One key aspect of this paper is the measure of evolu- 
tion in the distribution properties of galaxy overdensities 
from a continuous volume sampled with the same selec- 
tion function over a wide redshift baseline. As our volume 
sampled is still limited, errors on the analysis presented 



in this paper are dominated by cosmic variance. The tech- 
nique presented here will be applied to a larger sample as 
the VVDS observational program progresses. 
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Appendix A 

Here we describe the application of the Wiener filter- 
ing technique to deconvolve the noise signature from the 
VVDS density map. We de-noise data in Fourier space, 
noting, however, that an equival ent filtering can be di- 
rectly applied in rea l space (e.g., iRvbicki fc Presst Il992t 
IZaroubi et ailll995|) . 

Let us assume that the observed smoothed density 
field #o( x ), and the true underlying density field <5t(x), 
smoothed on the same scale, are related via 



<5o(x) = *r(x) + c(x), 



(40) 



where e(x) is the local contribution from shot noise (see 
eq. ©)■ The Wiener filtered density field, in Fourier space, 
is 



5 F (k) =jr(k)5 (k), 



^(k) = 



(4(k)> 



(#(k)) + (27T)3P e (k) 



(41) 
(42) 



where brackets denote statistical averages and where 
P e (k) = (27r)~ 3 (|e 2 (k)|) is the power spectrum of the 
noise. Assuming ergodic conditions for the noise, we can 
derive its power spectrum as P e (k) = (27r) _3 |e(k)| 2 . 
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The Power spectrum of the underlying theoretical den- 
sity distribution of galaxies, smoothed with the window 
F and taking into account the VVDS geometrical con- 
straints, can be derived from eq. (J5J. Specifically he the- 
oretical overdensity field smoothed on a certain scale R, 
which is sampled by an idealized survey with no selection 
functions, is 



(r — r p ) * F 



R 



1 



(43) 



If we assume that this density field is periodic on same 
volume V, having, for example, the same geometry of the 
VVDS survey, its Fourier transform is 



,ik-r j3 



(Ft 



(44) 



where the mean theoretical galaxy density p may be esti- 
mated averaging over sufficiently large volumes the VVDS 
data corrected for selection functions and sampling rate 
(see eq. (©)■ 

In eq. I|44|) , rii represent the occupation numbers of the 
infinitesimal cells d,Yi in which the VVDS volume can be 
partitioned (n^=0 or 1) and the sum is intended over all 
the cells of the survey volume. The Fourier transform of 
the smoothing window function F with which the discon- 
tinuous galaxy density field is regularized is, in the case 
of a Top-Hat spherical smoothing filter, 



F k = 3 



sin(kR) cos{kR) 
(kR) 3 (kR) 2 



(45) 



In order to obtain an estimate of the quantity (8 2 ^(k)} 
which enters the Wiener filter definition fea. l42f> . we com- 
pute the ensemble average of the squares of the modulus 
of the Fourier amplitudes (eq. I44|) and obtain 



iki\ 



where 



W(k) = ~ J e lk r d 3 r 



(46) 



The mean value of the occupation number is given by 
(Peebles 1980) 

(mrij) = p 2 d 3 r l d 3 r :j [l + f(r 4 - r rf )] 
(n,) 2 = (m) =pd 3 ri 

where the correlation function may be expressed via the 
Fourier conjugates 



ikr; 



(47) 



P(k) 



J d 3 r£(r)e ik r 



(48) 



We convert the sums into integrals over the occupation 
cells taking into account the specific VVDS geometry and 
obtain 



{S 2 T (k)) = \F k \ 2 Jd 3 k'P(k')G(k-k')^ 
+ \F 2 \(pV)- 1 + \W k \ 2 (l~F k )i 



were 



G(k - k') 



(2tt)= 



■\W(k-k 



'\|2 



(49) 
(50) 



(51) 



Note that, in the idealized conditions of a survey of 
infinite spatial extension (V — > 00) and no smoothing ap- 
plied to data the previous expression reduces to 



(SUk)) = (2ir) 3 P(k')5 D {k~k'). 



(52) 



Given the quasi pencil-beam nature of the first sea- 
son VVDS data, however, the sample of the 5 field can 
be described, with good approximation, as confined to a 
cylinder of length L (aligned along the redshift direction) 
and radius R. In this case the Fourier transform of the 
survey window function (cfr. eq. I46f) is given by 



fL, \ n Ji(k±R) 
W k = ]Q [-k ] {)2- n 



k,R 



(53) 



where jo and Ji are the spherical and first kind Bessel 
functions and where k — ^Jk 2 + k\ . 

The convolution integral on the right hand side of 
equation |5D1 can then be evaluated as follows (e.g., Kaiser 
& Peacock 1991, Fisher et al. 1993) 

J d 3 k'P(k')G(k-k') = JL J°° dyf {y-k^)F{y) 
where V=irR 2 L is the volume of the cylinder and where 



FW = f p ( 



\/x 2 + . 
R~ 



J 2 (x>) 



dx'. 



(54) 



Note that the theoretically expected variance of the 
density field in the VVDS volume is derived in real co- 
moving space. Thus, we have corrected the theoretical pre- 
dictions (cfr. eq. EOl in order to take into account redshift 
space distortions induced by galaxy peculiar velocities (see 
discussion in §5). 

We compute the galaxy density field on a regular 
Cartesian grid of spacing 0.5 /i _1 Mpc using the smooth- 
ing scheme presented in eq. JSJ). The resulting 3D density 
map is then Fourier transformed in redshift slices hav- 
ing line-of-sight dimensions dz=0.1. This partition strat- 
egy is implemented in order to describe consistently, using 
cylindrical approximations, the deep survey volume of the 
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Fig. 17. The overdensity distribution reconstructed on a 
scale R=8 /i _1 Mpc and in two different redshift intervals 
using the GALICS semi-analytical simulation (see §4.1) 
is plotted as a function of log(l + 8 g ). The dashed line 
represent the distribution of overdensities traced by sim- 
ulated galaxies with M. < — 15 + 51og/i. The shaded area 
(dotted line) represent the observed overdensity distribu- 
tion i.e., the distribution recovered after applying to the 
simulation the VVDS flux limit (I<24) and all the under- 
lying VVDS instrumental selection effects (see §4.1). The 
solid line represent the distribution of density contrasts 
after correcting the observed distribution with the Wiener 
technique. 

VVDS (whose comoving transversal dimensions are an in- 
creasing function of distance.) The Wiener filter at each 
wave- vector position is then computed by using eq. I|50(l . 
Finally, as described in §3, we select only the Wiener fil- 
tered density fluctuations recovered in spheres having at 
least 70% of their volume in the 4-passes, VVDS-02h-4 
field. 

In Fig. E|we use the GALICS semi-analytical simula- 
tion (see §4.1), to show the effect of the Wiener filter on the 
reconstructed galaxy density field. As Fig. El shows, the 
net effect of the correction is to shift towards low-values 
the density contrasts having low signal-to-noise ratio. By 
definition, in fact, F is always smaller than unity. Fig. 
IT7I shows that i) at the same density, the effects of the 
correction are bigger at high redshift where the density 
field is noisier due to the increasing sparseness of a flux- 
limited sample, and ii) at the same redshift the Wiener 
filter mostly affects the low-density tail of the distribu- 
tion where the counts within the TH window are small. It 
is also evident from Fig. 1171 that, in the density interval 
— 1 < log(l + S g ) < 1, the Wiener filtered distribution of- 
fers a better approximation of the underlying PDF, than 
the observed (uncorrected) overdensity distribution. 
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